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Abstract 

The main result of these notes is an analytical expression for 
the partition function of the circular brane model [1] for arbi- 
trary values of the topological angle. The model has impor- 
tant applications in condensed matter physics. It is related to 
the dissipative rotator (Ambegaokar-Eckern-Schon) model [2] 
and describes a "weakly blocked" quantum dot with an infi- 
nite number of tunneling channels under a finite gate voltage 
bias. A numerical check of the analytical solution by means 
of Monte Carlo simulations has been performed recently in 
Ref. [3]. To derive the main result we study the so-called 
boundary parafermionic sine-Gordon model. The latter is of 
certain interest to condensed matter applications, namely as 
a toy model for a point junction in the multichannel quantum 
wire [4]. 
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1 Introduction 



Since the seminal work by Fateev and Zamolodchikov [5] CFT models with parafermionic 
symmetry were extensively studied. In the original formulation the parafermionic current al- 
gebra appears as an extended conformal symmetry of self-dual multicritical points of sym- 
metric statistical systems [6-9]. Later this algebra has been employed in string theory [10]. 
Much effort has been devoted to the parafermionic models with a boundary [11-15]. 

In conformal field theory there is a special class of Conformally invariant Boundary 
Conditions (CBC's) [11]. In the case of parafermionic models some of CBC's are easily 
visualized in terms of the original symmetric statistical systems [11-13]. The fluctuating 
variables ( "spins" ) in such systems can be thought of as a set of k special points on the unit 
circle, 

cre{cj"}^~Ji with u = . (1.1) 

The simplest microscopic boundary conditions are those for which all spin boundary values 
are the same, say as = cu". The scaling limit of the self-dual multicritical system depicted 
in Fig. ^ is described in terms of the minimal parafermionic model [5] on the half line x < 
with the so-called fixed CBC imposed at x = 0. With some abuse of notation we refer to 
these CBC's, as well as the corresponding RG fixed points, as Bn^n = 0, 1, ... /c — 1). 



X 



Figure 1: Z^ symmetric statistical system with fixed boundary condition. The spins a fll.lj) 
are located at the vertices of the half-infinite square lattice. For an explicit form of self-dual 
Boltzmann weights see Ref. [9]. 



Yet another simple type of CBC is the "free" one. In this case, the microscopic boundary 
spins are free to take any value from the set (ll.lj) . In the present work, we denote such CBC 
and the associated RG fixed point by Sfree- The universal ratio [12-15], 



gfree 
gfixed 



Vk , (1.2) 



of the corresponding boundary degeneracies [16] is greater than one. Thus a unitary bound- 
ary RG flow TZ^'^^ from Bi^ee to Bn,n seems to exist, or at least does not contradict the 
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g-theorem [16,17]. Indeed, IZ^'^^ (the Ising model case) is a textbook example of boundary 
flow [18]. For an arbitrary k the existence of T^'-'^-* can be advocated within the general 
approach [20]. This boundary flow has been also argued in Refs. [14,19]. 

In this paper we study the QFT model underlying the boundary flow TZ^^\ In Sec- 
tion 12. II we introduce its Hamiltonian Hg . It is crucial that Hg depends on an additional 
dimensionless angular parameter, 6 = 6 + 27Tk. The Hamiltonian Hg''^ suffers from a spe- 
cific ultraviolet divergence, which is very similar to the "small instanton" divergence in the 
2D 0(3) nonlinear cr-model [21,22]. In Section [2.21 we consider a certain non-perturbative 
regularization of Hg''^ involving an additional Bose field which does not suffer from this prob- 
lem. The regularized model is referred to as the Boundary Parafermionic (BP) sinh-Gordon 
model. The sinh-Gordon parameter b plays the role of a regularization parameter and the di- 
vergence shows up when 6 — >■ 0. In Sect ion IT^ we propose an analytical expression ()2.36p for 
the partition function of the BP sinh-Gordon model in terms of solutions of a certain second 
order linear differential equation. Our motivation behind Eq. ()2.36p follows closely along the 
lines of Al.B. Zamolodchikov's unpublished notes [23]. Namely, in a view of Appendices A 
and B, the proposed formula is a straightforward generalization of Zamolodchikov's result for 
the conventional (without parafermions) boundary sinh-Gordon model. Unfortunately, even 
though Eq. ()2.36|) can be easily guessed, its derivation from first principles is still lacking. 
For this reason we perform some consistency checks for ()2.36|) in Appendix C. 

The main quantity we are interested in, is the partition function Zg'^^ corresponding to 
the Hamiltonian Hg''\ It can be extracted from the limiting 6^ — > +0 behavior of the BP 
sinh-Gordon partition function. An analytical expression for Zg''^ is given in Section 12.41 
Using this result, we describe a qualitative picture of the boundary flow governed by Hg^^ in 
Section ITHl It turns out that for each 6 from the open segment {2n — 1)tt < 6 < {2n + 1) vr, 
the boundary flow is terminated at the flxed point Bn,n- At the same time the RG trajectories 
associated with 6 = 7r(2n + 1) with n = 1, ... /c — 1 posses nontrivial infrared flxed points 
whose boundary degeneracies are all the same and equal to 2 cos(^^) gsxed- Hence, the 
infrared physics of the boundary flow TZ^'^^ essentially depends on the parameter 6. For this 

(k) 

reason we shall include an additional symbol in the boundary flow notation: TZg . 

In the rest of the main body of the paper we focus on the large-fc limit. In Section El we 
identify Zg = limfc^oo Zg''^ with the partition function of the circular brane model [1], and 
the parameter 6 is understood as the topological angle. In Section HI a particular ultraviolet 
regularization of the circular brane model is considered. The regularized theory is known as 
the dissipative quantum rotator model [2] and commonly used in describing the Coulomb 
charging in quantum dots. The applicability of our results to the quantum dot in weakly 
blockaded regime is briefly discussed in Conclusion. 

Numerous technical details of the work are relegated to appendices. In Appendices A 
and B we study the model which, from a very formal point of view, is the BP sinh-Gordon 
model considered at purely imaginary values of the sinh-Gordon parameter b. In essence, 
these appendices constitute work done in 1999 and which was never published before. It 
follows closely the approach developed in the series of papers [24-28]. In particular, some 
of the key statements involved have been already proven. In those cases we omit the proofs 
and refer the reader for details to [26-28]. Among results obtained in Appendices A and 
B is equation ()A.80|1 . which is interesting in its own right. It generalizes a conjecture of 



3 



p. Fendley and H. Saleur [4] for the DC conductance in the multichannel quantum wire 
to the case of finite temperature. It should be emphasized that, although the style of our 
presentation is somewhat sketchy, it seems possible to transform the derivation of ()A.80|) 
into a rigorous mathematical proof. 

In Appendix C we perform some consistency checks of the proposed exact analytical 
expression for the partition function of the BP sinh-Gordon model. Finally, for reference 
purposes, we collect some formulas concerning the low-temperature expansion of the partition 
function Zg'^^ in Appendix D. 

2 Boundary flow izf"^ 
2.1 Hamiltonian H^^^ 

Let us consider the minimal parafermionic model on a half lane, a; < 0, constrained by 
the free CBC By employing the magic of modular transformation [13] to the parafermionic 
characters [10,29], one can observe operators of the scaling dimensions 

A„ = ^^i^f^ {n = 0,...,k-l) (2.1) 

k 

among boundary fields of the theory. Below they are denoted as 

{^n{r)}tl . (2.2) 

The variable r labels points along the boundary a; = (see Fig. Notice that ^'o coincides 
with the unit operator /. We also reserve special notations for the two relevant boundary 
fields of the lowest nontrivial scaling dimension Ai = A^-i = 1 — ^, namely 

= ^ =^k-i. (2.3) 

The boundary fields ^^(t) can be thought of as Fateev-Zamolodchikov chiral parafermionic 
currents brought to the boundary point t} Thus, their boundary Operator Product Expan- 
sions (OPE's) have a form similar to those of the chiral parafermionic algebra: 

*,(ri) *„(r2) = I *,+™(r2) + . . . . (2.4) 

Here the sum j + m should be understood modulo k and the structure constants are the 
same as in [5]. With a properly normalized ^j, the structure constants read explicitly as 



_ I jj + m)\{k - j)\{k - m)\ 

j\m\k\ik-m~j)\ • ^ ' 

Let n\Z and H):' be the spaces of states of minimal parafermionic models on the half 
line constrained, respectively, by the free (i3free) and fixed (i3n,n) CBC's. Obviously, all the 
spaces Tin^ are formally isomorphic for different integers n: 

nl'^ ~ nP ~ . . . ~ n^l, . (2.6) 



precise meaning of this statement is discussed in Section IR21 
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By virtue of the simple microscopic nature of the fixed and free CBC's (see Introduction), 
it is expected that 

^tU©n=o^i'^ ■ (2.7) 

Then, the fields ^±{t), being considered as operators acting in ?i[^gg, intertwine the linear 

subspaces Hn^ fl'^free '^i±i(modfc) fl'^free- Heuristically, one can think about ^'±(r) as 
operators creating discontinuities in the fixed boundary condition Bn.,n — ^ Bn±i,n±i at the 
boundary point r (see heuristic Fig|2I). 
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Figure 2: Discontinuity diagram in the fixed boundary conditions corresponding to the 
insertion of the boundary parafermions. 

In this work we study the model described by the Hamiltonian 

Hi'^ = H[;i^f^[e'-l^^ + e^T^.] , (2.8) 

(k) 

where H^^^^ is the Hamiltonian of the minimal parafermionic model subject to the free CBC, 
*± = *±(0), and, hence, 

Hf^ : ni^l n^l . (2.9) 

The parameter in ()2.8|) carries the dimension [length]"^ , so that yu'^ sets the physical 
energy scale (the Kondo temperature) in the theory. It is useful to keep in mind that the 
OPE's (12 .41) . (j2.5|) specify the normalization of the boundary fields (j2.2|) modulo global 
transformations, 

^ uj"^ (a = 1, ... A; - 1) . (2.10) 

For this reason the angular parameter 6 = 9+2tt k in Eq. ()2.8|) is not completely unambiguous; 
it is defined modulo transformations 

6 ^9-271 a (a = l,...fc-l) . (2.11) 
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The partition function 



Tr, 



n 



(fc) 



exp ( - Hf^/T) 



(2.12) 



is an even 27r-periodic function of 6. Therefore the uncertainty in 6 does not affect the 
thermodynamics of the model, and we may consider 6 in the domain < 9 < n only. 



2.2 The BP sinh-Gordon model 

The key idea behind our calculation of ()2.12j) is that -ffg*^^ appears in a certain limit of a 
more general Hamiltonian, involving additional bosonic degrees of freedom. The generalized 
model is referred to below as the Boundary Parafermionic (BP) sinh-Gordon model. 
Let $(a;) and Il{x) be Bose field operators obeying canonical commutation relations 

[U{x),^{x')] = -2m6{x-x') . (2.13) 

The Hamiltonian of the BP sinh-Gordon model is described as follows: 

Hi%, = H\^l + H^^ + /. - ( e^^- + e-^*- ) , (2.14) 

where 

<i = i- j\x{U^ + {d.^r). (2.15) 

We do not impose any constraint on the boundary values of the canonical fields ()2.13p . Hence 
the Hamiltonian ()2.14p acts in the space 

^ = ^il®^i!^ (2.16) 

where H^^^^ is the same as in Eq. ()2.7|) and Hf*g is the space of states of the one component 
boson ^{x) with no constraint at the boundary x = 0. The parameter h in ()2.14j) plays 
the role of an external field coupled to the boundary values $b = ^{x)\x=q. To give an 
unambiguous definition of the coupling constant yU, one should specify the normalization of 
the boundary fields 

V± = *± e^^*^ . (2.17) 

In what follows these composite fields are assumed to be canonically normalized with respect 
their short- distance behavior, i.e. 

V±(ri) V^(r2) = \n-r2\-^''xl + ... , (2.18) 

and 

V^{n)V^{r,) . . . V^in) = ^ n - ^^-1'^'"'^ X + • • • ' (2-19) 
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while 

g±6v/fc*s(^^) g±fev^*s(^2) = _ ^^\2kh^ Xl+... . (2.20) 

Here 

1 62 

rf^l---- (2.21) 

is the scahng dimension of the boundary fields V± ()2.17|) . Also notice that the factor k\/k^^'^ 
in (j2.19|) is a product 11^=1 ^ij ^^e structure constants (|2.5p . 

Formally, when the dimensionless parameter 6^ in ()2.14|) is equal to zero, the Hamiltonian 
-^bshg becomes the sum of a free Bose Hamiltonian and H^q^ . In fact, because of the 
particular divergence generated by the interaction term in ()2.14j) at 6^ — > 0, the limiting 
procedure is not a trivial matter. The divergence appears due to the leading term in the 
OPE ()2.19j) . Keeping this in mind, we rewrite the Hamiltonian fl2.14|l as a sum of two 
operators 

Mlg = ^fi[*]+^2[*,$B] , (2.22) 

with 

H^m=Hfl + h^B + ^{Acos\i{h^^B)-Bl) (2.23) 

and 

/f2[*, ^b] = H^iL - (*+eVS*^ + *_e-Vl*^) - ^ (A cosh (6v^$b) - Bl) . (2.24) 

Our basic assumption is that the coefficients A and B for the local boundary counterterms 
in ()2.24|) can be adjusted in such a way that the limit of the operator H2 does exist when 
6^ +0. In other words, we expect that the singular behavior at 6^ —>■ +0 is fully controlled 
by the operator Hi ()2.23|) . which is just a Hamiltonian of the conventional boundary sinh- 
Gordon model. 

The 6^ +0 limit is the classical limit for Hi. Indeed, after the field redefinition cp = 
b\/k $ this becomes particularly striking. Since quantum fiuctuations of (p are suppressed as 
6^ — >• we may apply the saddle point approximation to account for their contribution. The 
bulk Hamiltonian for cp is free and massless, and, consequently, the saddle point is achieved 
for some constant classical field configuration which we denote as 16. Now, since all matrix 
elements of H2 are expected to be finite at 6^ = 0, we can safely replace in (j2.24j) by 
the c- number Obviously, the limit 

should be treated as the regularized version of the Hamiltonian H^q\ 

The saddle point configuration $ = corresponds to a minimum of the boundary 
potential C/[$b], which is a sum of the last two terms in the boundary sinh-Gordon Hamil- 
tonian ()2.23j) . The boundary potential is minimized at 

h 

sm{9) = i bVk — . (2.26) 
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We should also take into account the effect of Gaussian fluctuations around the classical 
saddle-point configuration. Let us write = i + Then, expanding the bound- 

ary potential near its minima, one gets 

U[^B] = ik {^S{9)-B) I + Acos{9) {6^B)' + 0{{6^Bf) , (2.27) 

with 

S{e) = 9 sin(^) + cos(e) . (2.28) 

Let 2'Gauss be the partition function of the Gaussian theory with the quadratic boundary 
potential, Z^\k) be the partition function corresponding to the regularized Hamiltonian 
(j2.25|) . and Z^^j^ (h) be the partition function of the BP sinh-Gordon model: 



ZZ,{h) = Trn [ exp ( - H\^ijT) J , (2.29) 

where the trace is taken over the space of states ()2.16|) . Then, the above consideration 
implies the following relation between these three partition functions as 6^ —>■ +0: 

--^)L^^+0 --"^ ^Causs . (2.30) 

Finally, according to Ref. [32], 

r(l + 2|cOS^ /Cx-2#cosW 
^Gauss = gD , 7^ • (2.31) 



47r|cos(^) 

Here gu = 2~3 is the boundary degeneracy associated with the Dirichlet CBC of the un- 
compactified Rose field and C is some nonuniversal dimensionful constant. 

2.3 Exact expression for the partition function Z^^l^^{h) 

Appendixes A-C are meant to demonstrate that the partition function ()2.29p can be ex- 
pressed in terms of the Wronskian Vr[6+, 6_ ] = 6+9u6_ — 6_(9nO+ of two particular 
solutions of the Schrodinger equation 

' -dl + ^^e-^ + e^)' Q{u)=0, (2.32) 

where 

Q = b + b-^. (2.33) 

More precisely, let 9+ and 9_ be regular solutions of (|2.32p as m +oo and u — > — oo, 
respectively. They are specified unambiguously via their own WKB asymptotics: 

0±(m) ^ (2/t)"5 exp (F(6^^ I ±m)) as m ^ ±oo , (2.34) 
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where F is expressed in terms of the conventional hypergeometric function, 



ku 

Q2bQ 2J 



P ( _ k k . 1 _ k I p~"^ 

4bQ k ^ ■ -i-Tiy 26Q' 2 ' 2bQ i ^ ) ' 

Then, the partition function ()2.29j) is given by^ 

The parameters of the Hamiltonian ()2.14|1 are related to ^ and k in ()2.H2j) as follows: 

^^k h 



2Q T 
1 k 



2'kT bQ 



2TTfl 



(2.35) 
(2.36) 

(2.37) 
(2.38) 



Notice that the Wronskian VF[0+ , 0_ ] can be viewed as a spectral determinant of the 
Schrodinger operator ()2.32|) . Indeed the spectrum of ()2.32p is simple, discrete and positive: 



spect = {en}Zi (e'>d>o if 



n > m) 



and, as follows from the Bohr-Sommerfeld quantization condition. 



nk 



n 



logn 



l + 0( 



log(logn) 
log n 



as 



n 



oo 



(2.39) 



(2.40) 



Therefore the spectral determinant can be defined through the equation 



ra=l 



(2.41) 



where Dq{k) is some overall ^-independent factor. At the same time, iy[G+, G_ ] is an 
entire function of the variable which shares the same set of zeros as (jCTIl .^ Using the 
standard WKB methods [33], one can show that log iy[ 0+ , 9_ ] = o(^^) as oo. Hence, 
as follows from the Liouville theorem, the ratio D /W does not depend on ^, and we can 
always adjust the overall normalization factor Dq in ()2.41|) to make it equal to one, i.e.. 



[U] 



(2.42) 



•^Ea. (|2.36|l is a generalization of the original proposal from [23] to the case fc > 1. For k ~ 1 Al.B. 
Zamolodchikov presented a set of convincing arguments in support of (|2.36() . In particular, the Wronskian 
was expressed in terms of the solution of the massless thermodynamic Bethe ansatz equation 
associated with the boundary sinh-Gordon model H2.23|l . 

^Let be a zero of M^[8+ , 8_]. In this case the solution 0+(u) and 8_(m) are linearly dependent. 
It just means that the Schrodinger equation possesses a solution which is both regular at x — > +oo and 



—oo, hence e Spcct. It is also obvious that T4^[8+ , 8-] 



= 0. 
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2.4 Exact expression for the partition function zjl^^ 

With the exact expression for Zbshg(^), we can now check the semi classical behavior ()2.3()|1 . 

Consider the spectral determinant D{k, k sin 6) in the limit 6^ +0. There is no problem 
with the limit for the differential equation ()2.32j) and for the solution G_|_ ()2.34|) . Namely, 



0+(m) exp 
where is a solution of 



(fcf^ + c+) « sin 



2/7r(fc-l) N 



(m) as 



+0 



defined by the asymptotic condition 

g-(u+C+)K 



^4- M 



exp 



K 



-((1 + z)^ 



V2k(1 + e")4 
The constant c+ in (j2.43p reads explicitly 



1) 



as u 



-oo 



(2.43) 



(2.44) 



(2.45) 



(2.46) 



where ip{^) 



dz 



logr(z) and 'jE is Euler's constant. Notice that the argument of the 



exponential function in ()2.43|1 vanishes for any odd k. 

The asymptotic condition ()2.34j) for B_ is singular as +0. For this reason the 

limiting behavior of this solution is a slightly delicate issue. It can be analyzed using the 
WKB approximation for B_: 



e 



(wkb) 



U 



exp < K 



bk 

e 2Q 



bk 



(2.47) 



where V{x) is given by 

p(u) = y {e ^ + e^) -sm^{9) with sin(^) = ^/k 



(2.48) 



Notice that the subtraction term e~^^ in the integrand comes from the asymptotic condi- 
tions ()2.34|) and ensures convergence of the integral in ()2.47|) . Consider now the argument 
{ . . . } of the exponential in ()2.47|) as 6^ — > +0. It is easy to see that { . . . } ^ F{u) — -g^ kS{9), 
where S{6) is given by ()2.28j) . while F{u) — > uk,cos{6) as m ^ — oo. Then 



„ — C-KCOSf 

6 2_ ^g^g-^ 



Q kb^ 



as 6^ ^ +0 



where S_ is a solution of the differential equation ()2.44j) such that 

H_(u) ^ 6^"+'=-)"^°^'' as u^-oo , 
and c_ is an arbitrary constant. It will be convenient to choose 

c_ = ^(|)+7s + 21og2. 
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(2.49) 



(2.50) 



(2.51) 



We emphasize that H_ can be defined by means of the asymptotic condition (|2.5(jp for 
< 6 < ^ only. For ^ < 6 < it, the solution S_ grows at large negative u and the asymptotics 
()2.45p does not define S_ unambiguously. Fortunately, the function S„/r(l + 2k cos6') is an 
entire function of the complex variable ( = cos{6) for u < 0. So the solution H_('u) can be 
introduced within f < 6* < vr through the analytic continuation with respect to the variable 
6 from the domain < ^ < ^. 

In a view of relations ()2.38p . the singular behavior of D{k,, K,sin6) as — +0 matches 
exactly the structure ()2.30|) provided 



A = Tk. 



B = Tk sin 



2^7r(fc-l)- 



(2.52) 



Furthermore, comparing the pre-exponent in ()2.30|) with a similar factor for D{k,k sin i 
we finally obtain the partition function for the regularized Hamiltonian ()2.25j) 



Zi'^K) 



gfixed 



2k cos ( 



r(l + 2/tcos^) 



(2.53) 



Two important comments are in order here. First, Z'^^^ is the function of the dimensionless 
variable /t, which is in fact the inverse temperature measured in the units of the RG invariant 
scale E^: 



Y 



As follows from Eq. ()2.38j) considered at 6^ 



27r 







27r/i 



rri 



(2.54) 



(2.55) 



The parameter /i should be understood now as the coupling of the Hamiltonian ()2.25|) . 

Second, the partition function for the regularized Hamiltonian ()2.25|1 is defined modulo 
an overall factor e'^-'^'^°'^^, where c_ is an arbitrary constant. By fixing the value of c_ via 
Eq. ()2.51|) we have enclosed a particular normalization condition for the regularized ground 
state energy 



^f) = -limjTlogZr] . 



(fc)- 



T-»0 



(2.56) 



Namely, if c_ is chosen as in (|2.5ip . then E^^"* satisfies the condition 



dff- 



:E. 



(k) 



\e=o 



. 



(2.57) 



In the case of odd k the regularized ground state energy satisfying ()2.57p is unambiguously 
defined. In particular, the dimensionful constant 



E, 



1 9=0 



(2.58) 



can be viewed as a universal physical energy scale in the model ()2.8p . Since the constant B 
(|2.52j) does not vanish for even values of k, the regularized Hamiltonian (|2.25|) contains an 
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additional counterterm of unit operator. In this case only the difference Eg^^ — E^^l^ turns 



out to be a universal quantity. For even k we choose Eg. ()2.58|1 as an extra normalization 
condition for the regularized ground state energy. 

Whereas ^ satisfying (E37I) and (E3H|l is a universal scaling function, the ground 
state energy Ef^ for the original Hamiltonian H^^^ ()2.8p does depend on details of the 
regularization procedure. More precisely, the above consideration suggests (see Eqs. (|2.24|) 
and (j2.25j) ) the following general form for E^^"^ 



^ = -E^ 



) cos(^) - Lf sin^ (Z^) + ^ , (2.59) 



(k) 

where L\ ^ are some dimensionless nonuniversal constants which are expected to absorb all 
divergences in the theory ()2.8j) .^ Hence, the partition function ()2.12|1 has the form 



Zg = exp 



cosie) - «:4'^ sin^ (^i^) Z^;\f^) , (2.60) 



where zII'\k) is given by Eq. ()2.53|) . For the regularization described above = + 0(1) 

9 (k) 

with b +0. In general, the divergent constants L\ 2 are expressed in terms of some UV 
cutoff equipped with the theory. 

(k\ 

One should note that some particular cases of Zg have been already discussed in the 
literature. E.g., the case k = 2 was studied in Ref. [18]. The differential equation ()2.44j) 
for k = 2 can be transformed to Kummer's equation, so that Zg has the especially simple 
form, 



(«:) = ^n^'^'.^le,, {^^f . (2.61) 
r(i + 2kcos"'(|)) 

Eq. ()2.61|) is in complete agreement with the result of [18]. It should be also mentioned 
that for sin (6*) = the partition function is described in terms of the differential equation 
originally introduced in [31] in the context of the so-called paperclip boundary flow. The 
difference between ^ = and 6 = tt appears in a choice of S_-solution. Namely, and 
g^j-Q ^Y^Q Bloch-wave solutions of ()2.44|) with sin(^) = [34]. We refer the reader to 
those papers for a wealth of data about the differential equation (j2.44p with sin(6') = 0. In 
particular, the thermodynamic ansatz equations describing the partition function Zg'^'^ (|2.12j) 
for 6 = and 6 = tt can be found in Ref. [34]. 



2.5 Qualitative pattern of the boundary flow IZ^ 

Our prime interest here is in the high- and low-temperature behaviors of Zg''\ With ()2.53p 

it is straight-forward to check that the boundary entropy Sg^^ = (l + T ^) log Zg''\ has the 
following ^-independent limit as T ^ 00: 

lim (T) = log(gfree) . (2.62) 

T — ^00 

very similar phenomenon was discussed in [30] for the minimal parafermionic models perturbed by 
the parafermionic currents in the bulk. 
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At the same time 



Mm Si'\T) = M^'-'^ as 0<9<n ^^.es) 
log(g2) as 9 = n 



Here the boundary degeneracy g2 is given by the formula 



sfixed 



Mm) 

sin ' 



(2.64) 



with s = 2. Whereas (j2.62j) is just a simple consequence of the fact that TZq starts from 
BfTee, Eq. (j2.63|) is much less trivial. It shows that within < < tt the infrared behavior 
of the theory is defined by the fixed CBC, but for 6 = n the RG flow is terminated at some 
nontrivial fixed point. 

The CBC associated with the boundary degeneracy g2 can be easily identified. Recall 
that the minimal parafermionic model possesses the Cardy states [13] with the additional 
conformal symmetry, namely symmetry introduced in Refs. [35-37].^ They are in one- 
to-one correspondence with the k{k + l)/2 number of primary fields of the model. Let 
us denote the CBC's corresponding to the invariant Cardy states by Bi^n with < n < 
I < k — 1. Their boundary degeneracies are given by Eq. ()2.64|) with s = l — n + In 
particular the boundary degeneracies associated with Bq^i, . . . , Bk-2,k-i and B^-ifl are all 
the same and given by the constant g2 = 2 cos(^^) ggxed- 

Now, in view of the global invariance ()2.10|) . ()2.11|) . it is easy to imagine the whole 
qualitative picture of the boundary RG flow IZ^q"^ for < < 2Tik. In the case = 4 it 
is depicted in Fig. El Among the flow lines shown in Fig. El there are trajectories adjacent 
to the big circle associated with zero temperature. In Appendix D we present an exact 
formula ()D.3|) for the normalized ground state energy (j2.56|) . (j2.57j) . It turns out that E^^^ 
is a singular function at = tt (see formula ()D.2|) ): 

Ef^ oc 1^ - 7r|T^ , (2.65) 

where 

4 = ^. (2.66) 

It is worth noting that cig is the scaling dimension of the so-called first energy boundary 
operator, e [5]. Hence, an effective Hamiltonian, describing the system in a vicinity of 
the nontrivial infrared fixed points Bn+i,n , contains this relevant boundary field along with 
irrelevant conformal and descendents of e and /. Notice that the coupling constant 
for £ vanishes at ^ = vr, 37r . . ., so the leading low-temperature corrections come from the 
first descendent of e of the scaling dimension 1 + The low-temperature expansion of 
(|2.53|) for the radial trajectories 6 = it, Stt . . . was studied in Ref. [34]. More details about 
the low-temperature expansions of ^ for e ^ vr are presented in Appendix D. 

^See Section IB.2I for the brief review of Wfc symmetry and invariant CBC's in the minimal 
parafermionic models. 

^In string language [15], Bi,n are the so-called "A-branes". The fixed CBC's Bn,n are DO-branes coinciding 
with the points on the unit circle. The CBC's Bin with I > n are interpreted as Z31-branes, the chords 
stretched between the points uj" and ojK The "lightest B-brane" from [15] is identified with the free CBC 

Bfree ■ 
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e = 7c 




e = o 



33 



e=77c 



free 



Figure 3: The boundary flow TZq in the 2D plane of coupling constants. Polar coordinates 
(p, 0) in this plane are some properly chosen "running" couplings, p = p{T,9), = (l)(T,9), 
such that liniT^oo piT, 9) = and pmax = liniT-^o p{T, 9). The parameter 9 is constant along 
each trajectory at T 7^ 0. 

3 Circular brane partition function 

3.1 zf^^ in the limit k ^ 00 

The regularized partition function (j2.53p possesses a finite limit as — 00 provided the 
physical energy scale E^, is kept fixed. Indeed, let us shift the variable u in the differential 
equation ^U^i : 



u = y- \og{k) . 
Then, the large-A; limit brings ()2.44|) to the form, 

-dl + K^ exp ( e^ ) - ^2 sm\9) 1 E{y) = . 
The solutions S+ and S_ remain finite for any fixed y when k ^ 00. Thus 

Z0{k) = lim Z^^^k) 

A;— »oo 



(3.1) 



(3.2) 



(3.3) 



can be calculated using the same Eq. (|2.53|) with , S_ ] being the Wronskian of two 

solutions of the differential equation ()3.2j) which are fixed by the asymptotic conditions 



ft cos 6 (,Qg g 



as y ^ —00 



(3.4) 
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and 

^ (2/s;)-5 exp - i e^^ - k Ei( i e?^ ) as y ^ +00 , (3.5) 

with Ei(2;) = £^ ^e~^. Again, the asymptotic condition ()3.4|) is apphcable for < 6* < | 
only. Within the domain ^ < 6 < ir the solution S_ should be defined through the analytic 
continuation with respect to the variable 6. 

Let us assume now that k is odd and consider the partition function (|2.6U|) at k ^ 00. 
One has 

lim =e"'=°^(^)^ Zein) , (3.6) 



k 

k = 'i, 5 



where L is the limiting value of the nonuniversal constant As follows from ()2.60p . the 
large-fc limit taken among the even k may cause the result to differ from (j3.6p by an extra 
multiplicative factor e'*'^°"''*. Obviously, this subtlety does not affect the physical content of 
the limiting theory and we shall ignore it below. 

3.2 Circular brane model 

Here we identify the QFT describing the large- A; limit of the model ()2.8|) . 

For finite k the minimal parafermionic model defined on the semi-infinite cylinder r = 
T + 1/T, X < 0, contains two holomorphic iI)^{t^x) = if^±{T — ix) and two antiholomorphic 
^±{t, x) = 'j/'_|_(r-|-ix) currents of conformal dimensions Ai = Ai = 1 — ^. When /c — > cxd, the 
unlocal parafermionic currents turn into chiral spin-1 currents. The latter can be bosonized 
in terms of a pair of fields satisfying free massless equations of motion in the bulk a; < 0: 

ddX = ddY , (3.7) 
with 9=1 {dr + i dx) and d = \ {d^ — id^). Namely, 

lim -j/jj. = ^^dY ^ ^ dX ±idY ^ 



Hence, the large-/c limit of ()2.8|1 is described by the theory with the bulk Euclidean action 
yibuik = - / dr / dx {dXdX + dYdY) . (3.9) 

^ Jo J~oo 

In view of the heuristic picture of the boundary interaction in ()2.8|1 (see FigEj) one may expect 
the boundary values X|a;=o = Xb and l^|x=o = Yb satisfy some 0(2)-invariant boundary 
condition. This boundary condition has been already identified for 6' = 0, vr in Refs. [1,34]. 
In both cases the limit ()3.6p coincides with the partition function of the circular brane model. 
The bulk action of this two-dimensional model is indeed given by ()3.9|) . while the boundary 
values of the Bose fields are subjected to a nonlinear constraint 

X| + Y| = 1 . (3.10) 

9o 
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Let us recall some facts about the circular brane model [1]. Due to the nonlinear boundary 
condition ()3.10|) . this theory needs renormalization. It has to be equipped with the ultraviolet 
(UV) cut-off, and consistent removal of the UV divergences requires that the bare coupling 
go be given a dependence on the cut-off energy scale A, according to the Renormalization 
Group (RG) flow equation 

A^ = -29l~Agl + ... . (3.11) 

The leading two terms of the /3-function written down in ()3.1H) were computed in [38] 
and [39], and indeed agree with the more general calculations in [40] and [41]. Eq. p.llj) 
implies that 

oo 

~ A g,' (l + J2''m9o), (3.12) 

m=l 

where the symbol ~ stands for asymptotic equality. The coefficients Cm in the asymptotic 
series fl3.12|) are not universal, i.e., they depend on the details of the regularization. 

It is important that the general definition of the circular brane model involves an addi- 
tional parameter, the topological angle 6. The configuration space for the two component 
Bose field X^^ = {X, Y) consists of sectors characterized by an integer w, the number of 
times the boundary value {Xb, Yb) winds around the circle (j3.1Up when one goes around 
the boundary at x = 0. The contributions from the topological sectors can be weighted with 
the factors e^"'^. Thus, in general 



E 



Z^'^^ (3.13) 



where Z^"^^ is the path integral 

ZM = I VXVY e-^'^-'kl^'^l , (3.14) 



evaluated over the fields from the sector w only. As it was pointed out in [42],'' the functional 
integral ()3.14j) has a particular nonperturbative divergence due to the small-size instantons 
[45,46], which cannot be absorbed into the renormalization of the coupling constant go. Due 
to this effect, the partition function ()3.13|) is expected to have the form ()3.6|) . where 

^ = 2 log (A) , (3.15) 

with A some ultraviolet cut-off regularizing the small instanton divergence, and k the inverse 
temperature measured in units of the physical scale ()3.12|) .^ This suggests to extend the 

''The circular brane model shows many shiiilarities with the 0(3) nonlinear cr-model (or n-field) [43,44], 
where a very similar effect of the small instanton divergence was known for a while [21,22]. 

^Notice that the instanton cut-off A is different from the perturbative cut-off A in the RG flow equation 
H3.11|l . i.e. their ratio A/A is some non-universal function of the bare coupling go- In Section^lwe find this 
ratio for some particular UV regularization of the circular brane model. 
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results of Refs. [1,34] and identify (|3.6|) with the partition function of the circular brane 
model Zq where Q plays the role of the topological angle: 

Ze = exp 2k cos(^) log {^) ) Ze{K) . (3.16) 
3.3 High-temperature expansion 

Here we discuss the high-temperature behavior of the partition function Zg. 

First and foremost we recall a simple qualitative picture described in Refs. [1,31]. As 
follows from Eq. ()3.11|) . the circular brane model is a asymptotically free theory - at short 
distances the effective size of the circle fj3.10|) becomes large. For small but nonzero values 
of the "bare" coupling constant Qq we can locally ignore the curvature effects and write 
the Bose field X^^{t,x) = + 5X^(r, x). It makes sense to split the fluctuational part, 
6X'^, into the components normal and tangent to the circle ()3.10p at the point Xq . The 
components must satisfy the Dirichlet and the Neumann BC, respectively. Therefore the 
"microscopic" boundary degeneracy for go <^ 1 does not depend on 6 and coincides with the 
product of the boundary degeneracies associated with the Dirichlet (gn) and the Neumann 
CBC (gAT = gci where C = 2n/y/g^ is the full length of the circle (j3.1(Jj) ). The above 
consideration supplemented by the common wisdom of the renormalization group suggests 
that as -E* ^ T ^ A, the partition function must develop the following leading high- 
temperature behavior: 



gfixed 

7W) 



l + 0{g{T)) . (3.17) 



In writing the above equation, we replace the bare coupling qq in the "microscopic" boundary 
degeneracy g%/ ^/go by the "running" coupling constant g(T). As follows from the two-loop 
/3-function ()3.11|) . g = g(T) can be defined by the equation: 

K = ^ = g''e~^ . (3.18) 
Also, the factor ggxed in (|3.17|) should be understood as 



gfixed — gD 



2 (3.19) 



The exact partition function ()3.16|) indeed possesses the leading short distance behav- 
ior (|3.17p . The systematic high-temperature expansion for ^ = was obtained in Ref. [1]. 
That result can be easily generalized to 6' 7^ 0, yielding an expansion of the form 

„ 1 / A \ 2k cos 6 
fixed / \ „-Kcos6» 



V9 



( \ \ ZK COS f 
^) g--^''^' Y.^^ 'r.{g.e) , (3.20) 



n=0 



where the coefficients Zn{g,0) are power series in the running coupling constant ()3.18jl . To 
describe some coefficients in fl3.20|l explicitly, it is convenient to represent Zg as 
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Here 



X = -T dl\og{Ze) 



and therefore the factor Zq is subject to the constraint 



89^ 



. 



(3.22) 



(3.23) 



It is possible to show that the topological susceptibility x expands explicitly as 



X^E,. [2 log (J) -\ogg + 3-fE + log 2 - 2-fEg + 0{g') - ^«:(l + 0{g)) + 0{k')\.{3M) 

The expansion coefficients of Zg ()3.21|1 are somewhat simpler than Zn in ()3.2()j) . For example, 

Zo = l-{l + ^E)9 + 0{g',K') , (3.25) 

and 

Zg/Zo = 1 - ^ (1 + 0{g')) sin^ (f ) + 0{k') . (3.26) 



Notice that the coefficients ^ and ^ in Eas. dT^ . (Unni) are in agreement with the result 
of two-instanton calculation from Ref. [47]. 



3.4 Low-temperature expansion 

For k = oo, the boundary degeneracy ggxed in ()2.fi3|l should be understood as in Eq. ()3.19j) . It 
implies that the circular brane boundary flow is terminated at the fixed point corresponding 
to the Dirichlet CBC = for any 6 within the domain < ^ < vr. At the same time, for 
6 = IT the flow possesses a nontrivial IR flxed point whose boundary degeneracy is given by 



g2 = 2 g; 



D 



(3.27) 



(see Eq. (|2.64|) with s = 2 and k = oo). The last equation suggests that the boundary values 
of are constrained to two points. We refer the reader to Ref. [34] for a comprehensive 
description of this nontrivial RG flxed point. 

Here we describe the low-temperature expansion for 6 ^ n. In this case the free energy, 
Fg = —TlogZg, admits an asymptotic expansion in terms of T^: 



Fe ~ cos{e) log (A) + Eg - T log(g^) - ^^i^) 



T 
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1=1 



2/-1 



(3.28) 



The regularized ground state energy Eg is a large-fc limit of the function ()D.3|) . It explicitly 
reads 



Eg/E, = Co- {-fE + log 2) cos^ + 



(3.29) 



da; 



TT cos^(^) 



r(a; + 
rfa;) 



5F2(x + i, i, 1; 2 



' 2 



COS 



cos=^(^) X 3F2(x + l, 1, 1; 2, II cos^O)) 
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Here 3F2 is the generalized hypergeometric function and 

/•I Y(x--) 

Cq = -v^ i dx = -1.44142 .... (3.30) 

Jo r(a;) 

The function ()3.29|) is plotted in Fig. HI Notice that it develops a non-analytical behavior at 
e = n: 



7r(7r — 6) / TV — 6 
log(7r -9)^ 



Ee/E.=2+ ::'^:-''i, +o{ , :-\^ ) . (3.31) 




e 



Figure 4: The regularized ground state energy ()3.29|1 in the circular brane model. 



The first two coefficients Fi in the expansion ()3.28j) explicitly read as follows: 

e 



Fii0) 



12 sine ' 



and 



1 



720 sin^ 9 Jo 



dt 



6 log^ 



sin^ 
sint 



7 log 



sin^ 
sint 



(3.32) 



(3.33) 



4 The dissipative quantum rotator model 

The circular brane model has useful interpretation in terms of Brownian dynamics of a 
quantum rotator. It was noticed a while ago in Ref. [48] that the free massless bulk dynam- 
ics fl3.9p is equivalent to the Caldeira-Leggett model of a quantum thermostat [49]. Upon 
fixing the boundary values X^(r) and integrating out the bulk part of the field X^{t,x), 
Eqs. dSini), (nni|l reduce to 



3 fixed 



-■Adiss[f?]--Ae[r;] 



(4.1) 
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where 



and 



A,av] = 7^ / dr/ dr' . \ ' ^( , (4.2) 
25-0 Jo Jo sm^(7rT(r - r')) 



■n 

Mv] = -7r- I drdrVir) . (4.3) 



2tt 

The variable ri[T) is the angular field defined through 



{Xb,Yb) = (cost/, sinr/) . (4.4) 



The path integral ()4.1|) requires some UV regularization. Among a large variety of regular- 
izations there is one of special interest for applications [2,50]: 

Abqk[v] = -^dissiv] + Mv] + ^c[v] , (4.5) 

where 



Ac[v] = j^ / drr/^ (4.6) 

4-f^C Jo 

The functional ()4.5|1 amounts to the Caldeira-Leggett action for the dissipative quantum 
rotator. In the weak coupling regime {go ^ 1) the additional term ()4.fjj) just provides an 
explicit UV cutoff of the dissipative action ()4.2|) with the cutoff energy A ~ Ec/go- Now, 
once the regularization scheme is chosen, all coefficient q in the asymptotic expansion ()3.12j) 
are determined unambiguously within the standard perturbation theory, in particular [1] 

ci = -! vr^ (4.7) 

The additional term ()4.6p not only makes the circular brane model perturbatively well- 
defined, but also regularizes the small instanton divergence. For example, the one-instanton 
contribution to the topological susceptibility (|3.22p in the Gaussian approximation reads as 
follows [42,51]: 



(pert)^ C ^ 



f da A da* e~-^c 

/ ■ 1 :^ + 0{9o 

Jaa'Kl 2771 I - aa* 



(4.8) 



Here A(j is the functional ()4.6p evaluated for the one instanton solution. The exponential 

_^(i) 

factor e c makes the integral over the instanton moduli \a\ < 1 finite. Indeed, using the 
explicit form of the one instanton solution "q^^^ [45,46], 

exp iir]^^\T)) = .1^^ , where z = e^^^'^ , (4.9) 
^ 1 — a* z 
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one finds 



Now we can integrate over the instanton moduli explicitly: 



(pert) ^ ^ ^ ^ 



2 log (^^) + 0(^7o, T/i5;c) J . (4.11) 



It is instructive to compare ()4.11|) with the high-temperature expansion ()3.22p . In essence, 
Eq. ()3.22|) requires that the bare coupling constant expansion should be of the form 

^(pert) [21og(A) -log^o + 37i^ + log2 + 0((7o,T/A)] . (4.12) 

Combining Eqs. (|4.1ip and (j4.12|) . (|3.12j) yields the following relations between the perturba- 
tive cut-off (A), the small instanton cut-off (A), and the energy scale Ec in (|4.6p : 

A = ^ , (4.13) 

and 

5 Conclusion: Charge fluctuations in quantum dot 

The dissipative rotator model ()4.5j) was introduced in Ref. [2] as an effective field theory 
describing tunneling of quasiparticles between superconductors. Later it has been involved in 
a surge of activity in the study of a low-capacitance metallic island (quantum dot), connected 
to an outside lead by a tunnel junction [39,42,46,47,51-53]. A remarkable feature of such 
a device (see FiglHl borrowed from Ref. [53]) is that the essentially many-body phenomena 
may be described via a single variable rj (conjugate to the number of excess charges on the 
island). As was shown in Refs. [2,46], the fermionic degrees of freedom can be formally 
integrated out yielding the Matsubara effective action A^a for rj. The Coulomb part of this 
action coincides with Aq + Ac fl4.3|) . ()4.6|) provided the energy scale Ec is understood now 
as the single electron charging energy, i.e., 

Ec=^, (5.1) 
where C is the island capacitance. The parameter 

% = — (5.2) 
Ztx 

has also the simple physical meaning of a dimensionless applied gate voltage V^: ng = 
CgVg/e, where Cg is the gate capacitance (see Fig. EJ. 



21 



2D EG 



dot 



Figure 5: Schematic view of a quantum dot connected to a bulk 2D electrode. The dot 
is formed by applying a negative voltage to the gates. Electrostatic conditions in the dot 
are controlled by the gate voltage Vg. Voltage Va applied to the auxiliary gates controls the 
transmission through the constriction. 



In addition to the Coulomb terms, the effective action contains a highly nontrivial dissi- 
pative piece which essentially depends on details of the island-lead interaction. In fact, the 
transmission coefficients { in the island-lead interface are coupling constants of Aes 

[N is a number of spinless electron channels). The most general effective action appears to 
be too complicated for practical purposes. It simplifies drastically in the special large- limit 
when each individual transmission coefficient Ta vanishes, but the dimensionless (measured 
in units e'^/h) tunneling conductance 



a 



1 

2^ 



N 



(5.3) 



a=l 



remains finite. In this limit ^lefr takes the form of the dissipative rotator action ()4.5|1 with 

1 



9o 



(5.4) 



so the basic thermodynamic quantities can be extracted from the corresponding partition 
function, Zdqr. The renormalized charging energy Ec 



(5.5) 



is of special interest. As it has been argued above, the quantum dissipative rotator model 
possesses the universal scaling behavior provided qq 1 and Ec ^ T. This corresponds 
to the regime of almost perfect transmission for the low-capacitance metallic island. In 
the scaling regime ^dqr turns into the circular brane partition function Zg. Thus the exact 
results obtained in this work give a complete description of the renormalized charging energy 
in the scaling limit. In particular at T = 0, 



E, 



c 



27r^E, 



T=0 



2 cos(e) log 



Ece- 



Att^E, 



+ 



^E 



a, E, 



c 



oo ) . 



(5.6) 



Recall that the function Eq / E* is given by Eq. ()3.29|) and plotted in Fig. |3J while 

E, = Eca^l-^ + 0(^-2) ) e"™' . (5.7) 
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It is worth noting that the quoted analytical results for Ec with ^ = have been numerically 
checked in Ref. [3] via Monte Carlo simulations of ^dqr- 

The quantitative description of the quantum dot with finite number of channels remains 
a challenging problem. For finite N, one would still expect some sort of universal scaling 
behavior characterized by the energy scale [46] 

JV 

E,(xEcl[ V^-Ta . (5.8) 

a=l 

In the case of = 2 this was explicitly demonstrated in Ref. [53]. Matveev argued that 
the universal behavior of the two-channel quantum dot is described by the Hamiltonian 
Hg^^ fj2.8p with k = 2. Thus the thermodynamic properties of the system can be extracted 
from the Chatergee-Zamolodchikov partition function ()2.6H) . With regard to = 2 and 
N = oo "integrable" cases, it would be interesting to explore the possibility of using H^''^ 
()2.8|) as a QFT Hamiltonian underlying the universal scaling behavior of the A^-channel 
quantum dot provided N = k. 
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A Appendix : The BP sine-Gordon model 

In this appendix we shall study the BP sine-Gordon model. It provides an essential step 
toward the main result ()2.53|) . 

A.l Definition of the model 

Let -fffree be a Hamiltonian of the boundary CFT model which contains two noninteract- 
ing sectors, namely the minimal parafermionic CFT and the Gaussian theory of the free 
uncompactified Bose field: 

if free = H^Jc + -f'^free • (^-1) 

Suppose the parameter b in Eq. ()2.2Hl is purely imaginary, 

b = if3, (A.2) 
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so that the boundary fields (|2.17|) take the form: 

V 

Their scaling dimension is given by 



V± = *± e^'^*^ . (A.3) 



1 /?2 

For < < 1, V± are relevant boundary fields normalized as in ()2.18|) . and in what follows 
we will focus on the model described by the Hamiltonian 

MS = ^free - f e^^* + e-^* V_ ) , (A.5) 



where and V are parameters. The parameter carries the dimension of [lengthY~^. 
The boundary interaction contains an explicit time dependence through the factors e^'^*. 
Below we always assume that > 0. Notice that, although by means of the formal change 
of the field variable, $ $ — t, the Hamiltonian ()A.5|) can be brought to the form 
()2.14p . effects of the external fields V in ()A.5|) and h in ()2.14|) are very different. Whereas 
the system fl2.14j) possesses the thermal equilibrium state even for non-zero values of /i, the 
system flA.5|) . at a nonzero V and a temperature T, develops a nonequilibrium steady state 
which can be thought of as the result of an infinite time evolution of the equilibrium state 
of the corresponding "free" system, with the interaction term adiabatically switched on. We 
denote by ( O )^^^^ the expectation value of an observable O over this nonequilibrium steady 
state. 



The model (jA.5j) was already discussed in the literature. In particular in the absence 
of the parafermionic sector {k = 1) it coincides with the massless boundary sine-Gordon 
model [55,56]. For this reason we shall call ()A.5j) the BP sine-Gordon model. Besides 
being an interesting model of Quantum Field Theory on its own, the theory finds important 
applications in condensed matter physics. In Ref. [4] the model ()A.5|) was suggested to 
describe a point-like impurity in the "multichannel quantum wire"^. In this context the 
integer k is the number of spinless electron channels, ao = k/S"^ is the conductance (in 
units e^/h) of the wire without impurity, while V is proportional to the voltage drop across 
the impurity. The quantities of special interest are correlation functions of the Heisenberg 
boundary operators 

V±(t)=e±^^*C/-i(t)V±C7(t), (A.6) 

where U{t) is the time evolution operator corresponding to the Hamiltonian ()A.5|) . For 
instance, the energy dissipation rate E for the steady state, produced by the time dependent 
boundary perturbation ()A.5|) . can be expressed in terms of the expectation values ( V± 
Indeed, as follows from ()A.5|) : 

E = -JoV with = i/i ( ( y+ )l% - { y_ ) . (A.7) 

Notice that the quantity 

aiV, T) = ao (1 + 27r f Jd/V ) (A.8) 

is interpreted as a DC conductance of the multichannel quantum wire [4]. 

^The cases k = 1 (without the parafermionic sector) and k = 2 were previously studied in Refs. [57-59] 
and [60,61], respectively. 
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A. 2 Perturbative expansions for {V±)\^^^ 

If > the system ()A.5|1 evolves towards a steady state which is characterized by the 
expectation values ( V± This is not a thermodynamic equilibrium state and no simple 
explicit expression for its density matrix is known. Nevertheless, one can describe the defini- 
tion of this state in terms of the real-time perturbation theory [54] . It is useful to represent 
the perturbative expansions for ( V± in the form 

{V^)i% = ^d,logA^, (A.9) 
where log A± are formal power series expansions, 

oo 

logA± = J] with s = /i^ (27rT)2'^-2 _ (A^^Q) 

n=l 

For k = 1 the compact formula for the perturbative coefficients a™'' were obtained in Refs. 
[27,62]. In essence the derivation from [27] can be applied without change to the BP sine- 
Gordon model. As a result, the following formula emerges: 



" ■ ■ n 



2n-l 

TT 



<:i,---,e2n-l j = l 



where the sum is taken over all arrangements of the "charges" ei, . . . , e2n--i = ±1 with zero 
total charge e + X]s=i ^ = 0, and 

i-i 

^. = e + 5^e, . (A.12) 

s=l 

The coefficients Je,ei...e2„-i in ()A.11|) are expressed through the (2n — l)-fold integrals over 
the real-time correlation functions: 

J,,,,,..,„, = / dti dt2--- dt2n-i e^^^A X (A.13) 



oo 

2n-l 



Here V^]^^^ (t) = e'*-^f'''=° V± Q-^iHivcc qj-q ^j^g boundary operators in the interaction represen- 
tation and 

((■■■»o-Tr^[---e-^]/Tr^[e-%-] (A.14) 

denotes an expectation value over the equilibrium thermal state of the free system. 
As follows from the real-time OPE 

Vi"*)(t)v|"*)(t') ^ (i(t-0)"''> 0<t-t'^0, (A.15) 
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the integrals ()A.13|) can be taken literally only if the scaling dimension < |. For ^ < d < 
1 we assume the "analytic regularization" common in conformal perturbation theory, i.e. 
the integrals should be understood as analytically continued from the convergence domain 
^ed < |. For the sake of illustration we consider explicitly the first coefficient a^^\ In this 
case the correlators appearing in the integrand of ()A.13|) are given by 



And 



V 



sinh(— vrTti 



2d 



(ti < 0) 



Calculating the elementary integral over ti, one obtains 



Att sm{7Td) 



r(l - 2d)T{d±2p) 
r(l -d±2p) 



with p = —i 



V 



(A.16) 



(A.17) 



A. 3 The system with g-oscillator 



Remarkably, the formal power series ()A.10|) can be interpreted as high - temperature ex- 
pansions of the equilibrium free energy for certain systems which differ from ()A.5jl in that 
they involve additional boundary degrees of freedom. Originally this was demonstrated in 
Ref. [27] for the boundary sine-Gordon model {k = 1), but the same reasoning can be apphed 
to the generalized model ()A.5|) . 
Let us define the Hamiltonian 



H 



free 



Vh 



/i 



(A.18) 



Here we use the same notations as in ()A.5j) . Formula ()A.18j) includes a new ingredient, 
namely, the operators h, a+, a_ which commute with all the parafermionic and bosonic 
degrees of freedom. They form among themselves the so-called "g-oscillator algebra" , 



[h,a±[ 



q a+a. 



q ^ a 



q-q 



-1 



(A.19) 



with 



Aird 



(A.20) 



Recall that the parameter 1 — ^<d<lis the scaling dimension of the boundary fields V± 
()A.4|) . Let p+ be some representation of ()A.19|) such that the spectrum of p+{h) is real and 
bounded from above. The Hamiltonian ()A.18|) acts in the space 



(A.21) 



with Ti. given by ()2.16|) . For > this Hamiltonian is bounded from below. Then, the sys- 
tem (fXT8|l possesses a thermal equilibrium state described by the standard density matrix. 

Strictly speaking the partition function 



e T 



(A.22) 
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is ill-defined even for the non-interacting system. Indeed, since the space of states has 
a structure of tensor product (see Eqs. ()2.16|) . ()A.21|) ). the partition function Z+{fi,V)\^=o 
factorizes into a product of three terms. The first one coincides with the boundary degeneracy 
gfree ()1-2|1 . The compoucut p+ in ()A.2H) gives rise to the factor Trp^[e~] which is well- 
defined if the spectrum of p+{h) is bounded from above. The problem comes from a divergent 
contribution of the zero-mode of the uncompactified Bose field subjected by the free CBC 
To cancel this formal divergence we shall consider below the ratio V)/Z+{0, V). 

Carrying out the finite temperature Matsubara procedure, one obtains the weak coupling 
expansion of the form 

f^TT^ = 1 + E ^ith . = /.^ {2nTr-^ , (A.23) 

where the expansion coefficients are given by 

4+)= r,...,,„ G,,...,,„ . (A.24) 

Whereas Te^.,,e2„ are expressed in terms of traces over the representation p+, 

_ Tr,Je^a_,,---a^,,J 

the coefficients Gei...e2„ are multiple integrals: 



T 



2n 



= / dn dr2---/ dr2„e^SA^.-. X (A.26) 
Jo Jo Jo 

Here ^^^i^^**^) (t) = 6'^"'^'''='= V-t e~'^''^'''^'='= are unperturbed Matsubara operators associated with 
the boundary operators ()A.3|) . and (( ■ ■ • ))o is defined in ()A.14|) . As well as in formula ()A.13|) . 
we should assume the analytic regularization of Gei...e2„ in the case ^ < d < 1. One might 
also note that all the coefficients Te^...e2j in ()A.24|) are determined unambiguously via the 
commutation relations ()A.19|) and the cyclic property of trace. Thus, they do not depend 
on a particular choice of the representation p+ provided Trp^ [e~ ] exist. 

Now, following along the lines of Ref. [27], it is possible to prove that the formal power 
series defined in equation ()A.10|) coincides with ()A.23|) . i.e., 

= . (A.27) 

Of course, A_ ()A.10|) is also related to a certain equilibrium-state partition function. Namely, 
let p_ be any representation of ()A.19|) such that the spectrum of p^{h) is bounded from 
below. Consider the Hamiltonian 



H_=H 



free 



+ Vh- p (^a+V+ + a_V.^ (A.28) 
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H_ 

Q- — 



(A.29) 



(A.30) 



acting in Ti.^ =71® p- and the associated partition function 

Z_(/i,\/) = Tr„_ 
Similarly to ()A.27j) . one can prove the relation 

- z_(o,v) • 

A. 4 Characteristic properties of A± 

Equation ()A.10|) shows that log A± are formal power series of the dimensionless variable 

s = ^? (27rT)2^-2 , (A.31) 

while the expansion coefficients ait"^ depend on the dimensionless ratio V/T. In what follows, 
we treat functions of the complex variable 



p = —\ 



V 

4^ 



(A.32) 



Here are the main properties of A±{s,p). 

• log A±{s,p) are formal power series of the form ()A.10|) . or, equivalently, 

oo 

A±(.,p) = l + 

n=l 



(A.33) 



The coefficients ait\p) ()A.10|1 and A^\p) ()A.33jl are meromorphic functions in the 
whole complex plane of p and 



AL'\p) = A^r^\-p) , a(~)(p)=aW(-p). 

Both (p) and A^t^ (p) are analytic in the right half-plane, 
^ep > Pn with 2 Pn = {1 — d) n — 1 . 
For n = 1, 2, . . . k, they have a simple pole located at p = p„ of residue 



Resp=p„[v4„(p)] = Resp=pja„(p)] 



(A.34) 



(A.35) 



(A.36) 



2tt 



.VkT{d) 



k\T{l-n{l- d)) 
2n\{k-n)\T{n{l-d)) ' 



There is no singularity at p = p„ for n > k, therefore the domain of analyticity ()A.35j) 
can be extended to the half plane p > pn — Sn with some finite 5„ > 0. 
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• The following leading power-low asymptotic holds in the domain of analyticity, 

a^n\p) Cn pl-2n(l-d) p ^ OO {^6 p > Pn) , (A.37) 

where C„ is some constant. 

• The formal power series A±{s,p) obey the so-called "quantum Wronskian" relation, 

e'^'''PA+{qs,p)A_{q-^s,p) - e~^'''PA+{q-^s,p)A.{qs,p) = 2i sin(27rp) , (A.38) 
with q = e^'^'^. 

The proof of functional relation ()A.38|) is based on results of the work [26]. It is outlined 
in Appendix B. All other properties readily follow from formulas ()A.11|) and ()A.23|) . For 
example the integrals ()A.26jl are entire functions of the variable p ()A.32jl . Thus, singularities 
of the coefficients A^\p) are due to the traces over the representation p+. A brief inspection 
of (fX^ shows that Al^\p) possesses only simple poles. They may be located only at points 
2p = —{d{n + l)+m) where / = 0, 1, . . . n — 1 and m = 0, ±1, ±2, . . .. Hence, aI^^ are 
meromorphic functions in the whole complex plane of p. The same, of course, holds true for 
the coefficients ai^\p). Furthermore, since ()A.34|) holds trivially for pure imaginary p, the 
claim is also true for any complex p. 

The analyticity of ai-t\p) in the right half-plane ^ep > pn for some p„ and the asymptotic 
formula ()A.37|) follow immediately from Eqs. ()A.11|) . ()A.13|) . To determine pn and calculate 
the residue one needs to examine these equations in more details. It is possible to show that 
the boundary of the domain of analyticity is defined by the integral J+_ei...e2„_i f|A.13|) whose 
"cumulative charges" rij ()A.12|1 have maximum admissible values^", i.e. 

€i = ... = tn-i = +1 , e„ = . . . = e2n-i = -1 • (A. 39) 

Moreover, as p ^ pn, the dominant contribution to this multiple integral comes from the 
region of integration where the operators V^"*^ with the same "charges" ej combine into 
two well-separated clusters: ~ ■■■ ~ |t„_i|T ~ 1 and |t„|T ~ ■■• ~ |t2n--i|^ ^ 1- 

Therefore, to determine the singular behavior as p — > p„ we invoke the cluster property and 
replace the correlation function in the integrand ()A.13|) by 

F{to, . . . ,t„-i) e™"(*"-*o) . . . ,t2n-i) , (A.40) 

where to = and^^ 

dn = n {n{d-l) + 1) (A.41) 

is the scaling dimension of the boundary operator ^'„e ^ ^ (see Fig.lH]). The latter is a 

^"Note that because of the factor n^"7^ sm{TTdrij) the sum (jA.lip enjoys a remarkable property, namely, 
it contains only the terms where none of the "cumulative charges" rjj (j = 1, . . . , 2n — 1) vanish. 
^^Here it is assumed that n < k. 
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4 




Figure 6: The amplitude defining the boundary p„ of the domain of analyticity "Rep > pn for 
a^\p) {n < k). The incoming and outgoing legs correspond to the insertions of boundary 
fields V^^^\t) and V^^^\t), respectively . For n < k the intermediate state carries the 
scaling dimension ()A.41|) and pn = 



product of the bosonic exponential field and the boundary field ^„ fl2.2|) . The correlation 
function appearing in ()A.40|) . 



F(to,...,tn-i) = (27rr)-'^"x 



(A.42) 



has a very simple form in the case n < k. Namely, for ti > . . . > t„ 

n-l 

F{h, ...,tn) = (27rT)"^ J] Ci„ J] ( 2i sinh {nTit^ - t,)) 



r(int) 



(n) 



-OO 



in=l 



2d- 2 



(A.43) 



m<j 



Here C 



(fe— m)(m+l) 



are the structure constants ()2.5|) with j = 1. Combining ()A.4njl 
with (|XTT|l one finds that the coefficients a„(p) (n = 1, . . . A;) possess a simple pole singularity 
at p = p„ = and 



Res„[aW(p)] =(-ir 



+1 



2^"-i7rfc!n! 
n'^k^{k — n)! 



n-l 



sin(7™rf)52 JJsin2(7rjrf) , (A.44) 



where S' is the Selberg integral [63, 64] 



S 



/Ml 
du2 
-OO 



n(n + l) 



TT 



r(rf) 



dM2n-i]^ 2 sinh 

T{l-n{l~d)) 
n{n - l)!n"ri sin(7rjd) 



Ui - Uj 



2d~2 



(A.45) 



Plugging ()A.45|1 into ()A.44|) one arrives at equation ()A.36|) . 

Finally it is easy to see that for n > k there is no pole at p = pn and a^\p) is analytic 
in some half plane Re p > pn — Sn where 5„ is a finite positive constant. 
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A. 5 Uniqueness of the solution of the Riemann-Hilbert problem 

The above-described properties of the formal power series expansions A±{s,p) constitute a 
certain Riemann-Hilbert problem defining these series uniquely. For k = 1 this was shown 
in the Appendix of the work [27]. Here we generalize the result of [27] to the case A; > 1. 

Let us substitute the formal power series (jA.lOjl into the quantum Wronskian relation 
()A.38jl . This yields an infinite set of relations for the meromorphic functions {p) of the 
following form 

sm{7ind + 2Txp) a^^\p) — sin(7m(i — 27rp) a^^\—p) = Rnip) , (A. 46) 

where Rn{p) are expressed through a[^'*(p) with / = 1, . . .n — 1 only. For example, 

Ri{p) = 0, (A.47) 
R2ip) = {qa[-'\p) + q'^a[+\-p) )' e^^'^ sin(27rj9)/2 . 
Introduce the function r„(p) such that 

then ()A.46|) takes the form 

rn{p) - rn{-p) = X (A.49) 

TT 

T{n-l-nd + 2p) T{n-l-nd- 2p) Rn{p) . 

In the case n = 1, the RHS in ()A.49|) vanishes and this equation supplemented by the 
analyticity condition ()A.35|) implies that ri{p) is an even entire function of the variable p. 
Moreover, in view of the asymptotic condition ()A.37|) . one should conclude that ri = const. 
The constant is defined by the residue condition ()A.36|) . This yields formula ()A.17|) . 

We shall prove the uniqueness of the solution of (|A.46|1 for n > 1 by using induction. Let 
us assume all the functions am\p) for m = 1, . . . n — 1 are already determined, so the LHS's 
of equations ()A.46|) and ()A.49|) are given. Suppose also for the moment that the complex 
parameter d belongs to the strip — 5 < ^ed < where 5 > is some small number. 
If 3fJe(^^ — d) is sufficiently small, then r„(p) is analytic in the half plane 3ftep > except 
for a simple pole at p = pn = ^ (^(1 ~ d) — 1). The residue of rn{p) at this pole is fixed 
by the condition ()A.36|) . Also, it is easy to see that r„(p) —> as p —>■ oo within SfJep > 0. 
Hence, relation ()A.49jl supplemented by the above quoted requirements for r„(p) constitutes 
a simple factorization problem. The uniqueness of this Riemann-Hilbert problem follows 
immediately from the Liouville theorem. To complete the induction step we note that the 
function r„(p) within ^ed < 1 can be obtained via the analytic continuation in the variable 
d from the strip — - S K'Sied < —. 

Thus we see that flA.46|) provides a recursion for the evaluation of ait\p) which allows 
one to express aii'\p) in terms of an {n — l)-fold integral. Explicit formulas for a^^^ and 
af^ in the case k = 1 can be found in the Appendix of [27]. Those formulas can be easily 
generalized to A; > 1. It is worth noting that such integral representations prove to be 
convenient for numerical calculations because the integrals converge very fast at infinity. 
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A. 6 Exact expressions for A± 

It was observed some years ago in the work [65], that A±{s,p) for k = 1 and some particular 
values of p can be related exactly to the eigenvalue problem of a certain ordinary differential 
operator. This remarkable observation was generalized and proven for any p in Ref. [28]. 
Here we extend the result of [28] to all integers k > 1. 
Let us consider the Schrodinger equation 

^^dl + ^[exp{-^)+ exp (^) f-e} e(w) = . (A.50) 

For < < 1 the potential term decays at -u ^ — oo and therefore for real C > there is a 
Jost solution which is asymptotically plane wave as m ^ — oo: 

e_(M,0 ^e-'«" . (A.51) 

As a matter of fact, this particular solution is a meromorphic function in the whole complex 
plane of the parameter ^, so it solves ()A.50|) for all values ^ except points where 6_(m,^) 
has simple poles (see, e.g., [66]). Notice that B_(m, — ^) is another, linearly independent (for 

^ ^ 0) solution of (IX3n|l . 

Since the potential in ()A.50|) grows rapidly as n ^ +oo, the Schrodinger equation admits 
also a solution which decays at large positive u. We denote this solution by 0+(u, ^) and fix 
its normalization by the condition 

e+(M,0 ^ (2/t)~5 e^^'^l") as u^+oo, (A.52) 

where F{b \ u) is defined in Eq. fj2.35p . 6+('U,^) is an entire function of and solves differ- 
ential equation ()A.50|) for all values of this complex parameter. 

The main objects of our interest are properly normalized Wronskians of the above- 
introduced solutions, namely, 

D±{n, = f(I^ 2ie(i~-^rf)) ^ ""^^ ~ d))^''^^'-"^ iy[e+(n,0,e_(«,±0] ■ (A.53) 
Here d is given by ()A.4|) and the overall factor provides the normalization condition 

limL)±(/t,0 = 1 • (A.54) 

K— ♦O 

In what follows we will show that for < < 1 and all values of n and ^ 

A±{s,p) = D±{K,i) , (A.55) 
provided the variables are identified according to the relations: 



d 1 — d 



27T 



VkT{d) 



(A.56) 



More specifically, we will check that D± considered as functions of the variables s and p, 
obey the same set of conditions ()A.33|) - ()A.37|) as A±{s,p). Thus ()A.55|) will follow from the 
uniqueness of the solution of the Riemann-Hilbert problem. 
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In order to examine properties of D± it is useful to make a change of the variable in 
(IX50ll . 

z = j-^ + \og{2Kil-d)) , (A.57) 

and rewrite the differential equation using the notations s, p ()A.56|) and d ()A.4|) : 

{-dl+ e^' + V + ^U{z) ) Q{z) = , (A.58) 

with 

7 1 / r)2— d \ 2m 

The solutions 9± ()A.5Hl . ()A.52jl should be considered now as functions of the complex 
parameters s, p provided 1 — < d < 1 and z is real. In fact, it is convenient for the 
present discussion to change slightly their normalizations: 

e4z,p) = {2Kii-d)f''^^'-''^ e_KO, 

e+{z,p) = -1= e+KO- (A.60) 

V 1 — « 

Then formula ()A.53|) takes the form: 

Di{s,p) = ^- Y^^^^^ W[Q+{z,p),Q_{z,±p)]. (A.61) 

Both G+ and 0_ are entire functions of s. In addition, 0_ is a meromorphic function 
of p, while 6+ is an even entire function of this complex parameter. Hence D±{s,p) 
can be expanded in a power series of s similar to ()A.33|) . and the expansion coefficients, 
Dlt\p) : d'^\p) = D^\—p), are meromorphic functions of p. 

To proceed further we will need the following conventional formula for D^{s,p), valid 
over the domain ^ep > 0, 

/oo 
dz hpie') 5U{z) x{z) . (A.62) 
-oo 



D+{s,p) 

Here Ip{x) is the modified Bessel function and xi^) solves the Lipman-Schwinger equa- 
tion, 

/oo 
dz' G{z,z')5U{z')x{z ) , (A.63) 
-oo 

where G{z^ z') is the Green function of ()A.58j) with 5U = subject of the asymptotic 
condition \imz-,oo G{z, z') = 0. Notice that x(z) differs from the solution 9+(z,p) 
by some overall ^-independent factor only. The coefficients dIt^^ (p) can be calculated 
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within standard perturbation theory developed with respect to 6U. In fact, to deter- 
mine dI^\p) for given n, one needs to perform n — 1-perturbative iterations in the 
Lipman-Schwinger equation. In particular, 

/oo 
dz hpie') K2p{e') 5U{z) +0{{5Uf) . (A.64) 
■oo 

Making use of ()A.59|) and the table integral 



oo 



one finds for n < k: 



D^+\p) = [ ^=—-] r(-l + n(l-d))x (A.66) 



k\ T{l-n{l-d))T{l-n{l-d) + 2p) 



1/ 



+ 



n\{k-n)\ 4:^T{n{l -d) + 2p) 

Here the dots mean the terms corresponding to the higher-order perturbative in 5V 
contributions. Further, one may argue that the first-order contribution (jA.66j) is an- 
alytic in the half-plane Qm{p) > Pn = ^{^{1 — d) — 1), whereas the omitted terms 
have a wider domain of analyticity. Hence, D^\p) with < A; is an analytic function 
for '^m{p) > Pn and has a simple pole located aX p = pn with the same residue as in 
()A.36j) . Similarly, it is possible to show that d'^\p) with n > k is analytic in the 
half-plane ^m{p) > pn — 6n with some 5„ > 0. 

The large p behavior of the expansion coefficients dl^\p) for the series logDli^\s,p) 
can be explored by means of standard semiclassical methods. A simple WKB analysis 
of the Schrodinger equation ()A.58j) shows that d^\p) — > C„ pi-2n(i-o!) ^ ^ g^,^^ 
Qmp > Pn, with the constants C„ given by 



C„, = (-1)"+^ k 



ZTT 



X (A.67) 



Vkr{d). 

r{n-{l-d)kn)T{l - {1 - d) n)T{-l + {1 - d) n) 
20Fn! r(l -{l-d)kn) ' 

Finally, let us prove that D±{s,p) ()A.61|) obey the quantum Wronskian relation. We 
start with an observation that the following transformations of the variable z and the 
parameters s and p, 

A: z z , s ^ s , p ^ —p , 

: z ^ z + in , s ^ g^^ s , p p , (A. 68) 

leave the differential equation ()A.58jl unchanged while acting nontrivially on its solu- 
tions. It has been mentioned already, that the transformation A applied to the solution 
Q_{z,p) ()A.60j) yields another solution, and the pair 

e_ = e_(2,p), Ae_ = e_(z,-p) , (a.69) 
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forms a basis in the space of solutions of ()A.58|) . Indeed, 

iy[e_,Ae_] = -4p , (a.to) 

therefore the solutions ()A.69|) are linearly independent provided p ^ 0. The solution 
0+ ()A.60|1 can be always expanded in the basis ()A.69|) . 

e+ = ci 0_ + C2 Ae_ . (A.71) 

Making use of ()A.70|) we have 

C2 = ^W{s,p), (A.72) 
with W{s,p) = W[Q+, 6_] . The solution 9+ is an even function of p, therefore 

c, = -j-Wis,-p) . (A.73) 
Let us apply now the transformation (l ()A.68j) to both sides of ()A.7Hl . Since 

ne4z,±p) = e±2^'pe_(2,±p), so 

ne+ = ^ ( e-^'^'P Wisq-^ p) Ae_ - e^^'^ W{sq~\ -p) e_) . (A.74) 

Now, it is not difficult to check that 

W,[Q+,nQ+] = -2i . (A.75) 
Combining formulae ()A.71|) . ()A.74|) and ()A.75|) one obtains the relation 

Q^^^PW{s,p)W{q~'^s,-p) -e-'^^'"PW{q-'^s,p)W{s,-p) = Sip . (A.76) 

Hence D±{s,p) ()A.61|) satisfy the quantum Wronskian relation ()A.38|) with A±{s,p) 
replaced by D±{s,p). 

Notice that once the relation ()A.55|) is established we may assert that series ()A.23|) . ()A.33|) 
have an infinite radius of convergence. 

A. 7 Formula for the DC conductance 

Since the parafermionic boundary sine-Gordon model ()A.5|1 has been of a certain interest in 
the quantum impurity problem [4] , we present here an exact formula for the DC conductance 
l8|) which follows immediately from the above consideration. 
Using equation ()A.9|) . the conductance can be written in the form 

a{V, T)/ao = ^ + ^ ^Aog . (A.77) 
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This formula has a nice interpretation in terms of the Schrodinger problem ()A.50|) . Indeed, 
as f — > —oo, the solution 9+(f,^) ()A.52|) develops the asymptotic behavior: 

Q+{u, ^ const ( e'«" + S{^, k) e''^" ) , (A.78) 

where 

is the reflection scattering amplitude for the Schrodinger equation ()A.50|) . Therefore, in view 
of ()A.53|) and ()A.55|) . the DC conductance is expressed in terms of S{^, n) as follows, 

a{V,T)/ao = l + ^j^^^^-^ ^d, log S{^, k) , V/T = ^ {1 - (3') . (A.80) 

Recall that = a^/k is the dimensionless conductance per channel of the wire without 
impurity, while k in ()A.80|) is the dimensionless inverse temperature measured in unit of the 

fc — QQ 

Kondo scale: k = E^:/T, oc /i . 

Formula ()A.80jl extends the result of [28] to the case k > 1. We note also that as T — 0, 
the dimensionless variable p = —iV/{4:7rT) oo. In this limit, the coefficients a^\p) 
develop the simple asymptotic behavior ()A.37|) where the C„ are given by (ESlI)- This 
immediately yields the large- \^ expansion of the conductance at T = proposed in Ref. [4]: 



a{V,0) ^ ^,„r(n + l-fcn7)r(l-n7) [47rE,7T'"^ 



ao 2 ' n!r(l - fcn7)r(| -n7) 

where 7 = 



V 



(A.81) 



B Appendix : Quantum Wronskian relation 

The purpose of this appendix is twofold. First, it is intended to sketch how the quan- 
tum Wronskian relation ()A.38|) can be proven for the perturbative series A± ()A.27|) . ()A.30|) . 
Second, it provides a short review of some well-known issues of CFT which are useful for 
understanding the models under consideration. 

Our previous discussion has always dealt with the so-called open string picture (channel). 
In this Hamiltonian picture the variable a; < plays the role of the space coordinate, while 
the coordinate r is treated as the compact Euclidean time. We now want to switch over to 
the "closed string" channel. In this channel the world-sheet coordinate x is viewed as the 
Euclidean time, so that the CFT lives in the finite volume 1/T. The closed string channel 
Hilbert space is embedded in a tensor product of the chiral spaces of states of right and 
left movers. Below the following notational convention is applied: If the open string Hilbert 
space was denoted by Ti, then the right and left movers' spaces of states in the closed string 
channel will be denoted (with some abuse of notations) as H and H respectively. Of course, 

the spaces H and 7i are isomorphic to each other. We will focus mainly on the "right" 
component H. 
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In the closed string channel an effect of the boundary at x = can be described in 
terms of the boundary state which incorporates all information about boundary conditions 
[13,55,67,68]. In the case of the "free" CFT ()A.1|) the corresponding conformal boundary 
state I B )free takcs the form. 



where 



i?)f..ee=|5)!fi®|5)il, (B.l) 



5);f2c7i:(*)®7i:(*) (B.2) 



is the boundary state of the Gaussian CFT subject to the von Neumann boundary condition, 
while 

|5);2^c7i('-)®7Y(^) (B.3) 

is the boundary state of the minimal parafermionic model with free boundary condition. 

As an important step toward the derivation of ()A.38|) . we will show that the power 
series A±{s,p) are vacuum eigenvalues of certain operators acting in the chiral Hilbert space 
7Y(*) ®7YW. With this aim in mind we need to describe explicitly W'^^ and 'H^^\ as well as 
the boundary states ()B.2j) and ()B.3|1 . 



B.l Chiral space 7^^*) and boundary state |5)[^g 

A general solution to the bulk equation of motion A$ = has the form: 

#(r, x) = 0(r - ix) + 0(r + ix) , (B.4) 

where (f){v) and (l>{v) are right and left chiral fields respectively. In the case of the uncom- 
pactified Bose field $(r + 1/T,x) = $(r, x), the chiral fields admit the Fourier mode 
expansions of the form: 



(f,{v) = i $0 + vrt;r Ho + i y ^ ^-2nimvT ^y = r-ix) , (B.5) 

^ — ^ m 



and 



Mv) = i $0 _ vrt;T Ho + i V — e'"'"*"^ (v = T + ix) . (B.6) 

m 

The canonical quantization procedure leads to the following nontrivial commutation rela- 
tions: 

[*o,no]=i; 0„J = 0„] = f (5„+„,o • (B.7) 

Let jFp (jFp) be the Fock space, i.e. the space generated by a free action of the operators 
0n (0n) with n < on the "vacuum" vector | P ) which satisfies 

0JP) = 0JP) = O, for n>0; 



no|P) = P|P), no|P)=P|P). (B.8) 
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Then the closed string space of states coincides with the direct integral f^J^p ^J^p, and the 
corresponding chiral Hilbert space is given by 

= f J^p . (B.9) 



The boundary state ()B.2|) associated with the von Neumann boundary condition obeys 
the equation, 

{dMr.x))l^^\B)^;i = Q. (B.IO) 
It reads explicitly as follows [67, 68] 

oo 

\B)tl = gD exp(5^|(-ir 0_„0_„)|vac)(*) , (B.ll) 

n=l 

where | vac)^*-* = | ) ® | ) is the ground state of the Gaussian theory in the closed string 
channel. Now, in view of ()B.11|) it is easy to see that for l/T>ri>--->r„>0 the 
following relation, involving the correlators of the boundary exponential fields, holds: 

i!i(5|e^'^^*nri) ■■■e''^"*^(T:„)|vac)(*)= (B.12) 
g-2™PE;=i-. ^ p I e2i'^i'^(ri) ■ ■ ■ e2^''"'^(r„) \P). 



Here 



g2ia0^^^ = (2mT)'^ exp (2a ^ e^"^™"^) exp (la $o + ^marT Hq) x 

m>0 ^ 

exp ( - 2a ^ ^ e-2.in^rTj ^B_^3^ 

m>0 ^ 

is an intertwining operator for the Fock spaces: 

e2-<^(r) : ^ ^p+, . (B.14) 

B.2 Chiral space ii.^^^ and boundary state 



free 



B.2.1 Parafermionic and current algebras 



Here we review shortly some well-known facts about the chiral Hilbert space of the minimal 
parafermionic models. The details can be found in Refs. [5,36,37]. 

Again, let the variable f be a complex coordinate, f = r — ia;, on the 2D cylinder. The 
chiral parafermionic algebra contains a set of currents {V'n(^)}„-i of conformal dimensions 
(|2.1|) . This algebra possesses a global symmetry with respect to the transformations 
ij^^ — * uj'^"' with a = 0, ... A; — 1, and its invariant component can be generated by 
the currents with the lowest conformal dimension, i.e., ij)^ = ij^i and = V'fc-i) through 
the OPE: 

^ (m - vf ( W^{u) + W^{v) )+...}. (B.15) 
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The spin-2 current W2 flB.15|) . in turn, generates the Virasoro algebra with the central 
charge 

The higher-order terms involve higher-spin currents, W3, W4. . . forming the algebra 
introduced in [36, 37] with the special value ()B.16jl of the central charge. Notice that the 
higher currents can also be produced recursively from singular parts of OPE's of the lower 
currents. In this sense the chiral algebra is generated by two basic currents W2 and 
W3. 

The parafermionic algebra has a set of irreducible representations (irreps), 

{Mj\ J = 0, i, ... 1} , (B.17) 

obtained through the action of the '0+ and V- modes on the highest vectors (see Ref. [5] 
for details). The irrep M. j decomposes into a direct sum of subspaces specified by the Zk 
charge i (number of ■?/'+ modes minus number of t/?. modes modulo k): 

Mj = ©,to Mf . (B.18) 

Each component Aif in ()B.18|1 is an irrep of the Wk algebra. 

Recall that an irrep of the algebra is obtained by factorizing a highest weight modules 
over submodules of "null- vectors" [36,37]. In the case under consideration, the highest weight 
can be thought of as a pair of eigenvalues (A, w) of mutually commuting operators 

1 jT 1 jT 

^ + (2vrT)-^ / ^ W^iv) , (2vrT)-2 / ^ W^iv) , (B.19) 
24 Jo Jo 

corresponding to the highest weight vector. In fact, it is most convenient to label the highest 
weight vectors by means of another pair of numbers (J, M) such that 



A 



(J,M) 



J(J + 1) M' 



2 



k + 2 k ' 

w^j,M) = ^{^-^^^M'-6J{J+l) + k) . (B.20) 

If W(j,M) denotes the irrep of highest weight (A(j^a./), ""^(j.m)) ()B.20j) . then the structure 
of the parafermionic irreps ()B.18j) is described as follows: 

mS'^ ^ W(o,o), ^i^'^ ^ {j = l,...k-l) , (B.21) 

M':^ - W(j,j_,) for i<2J (J > 0) , 

and 

Mf ~ mV , foT i-i' = 2J (mod k) . (B.22) 

2 
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Let S be a set containing ^^^^"^-^ number of pairs (J, M) of the form: 

E = {(J,M) I J = 0, |, ...|; M = -J,-J + l, ... , J; 

(J, -J) ^ (I - J, I _ J) } . (B.23) 

As implied above, the chiral parafermionic algebra can be thought of as an algebra of inter- 
twiners in the set 

{>VWM)}(,,^),, (B.24) 
of "integrable" irreps of W^. In particular, 

i^±{v) : W(j,M) ^ W(j,M±i) • (B.25) 

Notice that the pairs (J, — J) and (| — J, | — J) are subject to the equivalence relation in 
()B.23|) since the corresponding highest weights (A, w) ()B.20|) are the same, i.e., W(j,_j) ~ 
W(fe_jfe_j)- Finally, the closed string Hilbert space in the case of minimal parafermionic 

models is given by the direct sum ©(j,M)es VV(j,m) (S> W(j,m), while the corresponding chiral 
Hilbert space can be described as follows: 

n^'^ = ©(J,M)es W(j,M) . (B.26) 

B.2.2 Wfc invariant boundary states 

The Wfc algebra posses an automorphism of the form 

C: Ws^i-iyWs, s = 2, 3, .... (B.27) 

Thus [13,68], there are two types of invariant CBC's. In the string theory they are 
referred to as A and B-branes [15]. The corresponding boundary states obey the equations 

A-brane : [W s{r) - {-ly W s{r) ] B ) = , (B.28) 

and 

B-brane : [Ws{r) - Ws{t)]^^J B) = . (B.29) 

Let \2(j,m))a be the A-type Ishibashi state [13,15,68] corresponding to the integrable 
irreps W(j,m) from the set ()B.24|) . In other words, it is a solution of ()B.28|) in the space 
y^iJ,M) ® W(J,Af), that is unique, provided 



{ ( ( J, M) I ® ( ( J, M) I } • I J(j,M) )a = 1 , (B.30) 

where | {J,M)) is the highest weight vector of W(j^m)- Now, following the formalism of 
Ref. [13], we can construct a set of states satisfying both the Cardy consistency condition 
and dESH): 

giJ,M) 

1%^))= (B.31) 

(J,M)eE y 5(0,0) 
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Here S^^j'^^ is a modular matrix of the characters of the integrable irreps.^^ The common 

wisdom then is that the ^'•^^'''^ amount of Cardy states ()B.31|) are the boundary states 
corresponding to some local, invariant boundary conditions. In Section l275l these CBC's 
have been referred to as Bi^n (0 < n < Z < — 1), where the integers /, n are related to the 
pair (j, m) in ()B.31jl as follows: I = k — j — m and n = j ~ m. In particular, the k amount 
of the boundary states { | B^^k k_^-^ ) correspond to the k possible types of the fixed 

CBC's, while | | B^k^ >^-s) ) J s-o boundary states associated with the non-trivial 

fixed points of the boundary flow Tlg^\ 

Let us consider now equation ()B.29j) . Again, to construct nontrivial boundary states 
of the B-type one must draw on the B-type Ishibashi states \T{j^m))b, i-e., the unique 
normalized solutions of ()B.29|) in the space W(j,m) ® VV(j_m)- Obviously \T(^j^m))b exists 
only if J and M satisfy the condition 

Wij,M) = , (B.32) 

where the polynomial W(^j^m) is given by ()B.20j) . The following solutions of the Cardy con- 
sistency condition [13] can be easily verified [12-15]:^^ 

\B,) = VkY, -^|JW0))b (j = 0,1,2,...[|]) . (B.33) 

j=o Y'S'(o'o) 

There are strong reasons to believe [12-14] that the solution | Bj ) with j = corresponds 
to the free CBC: 

|i?)Se = |5o) . (B.34) 

By virtue of the above described structure of this boundary state, one can show that the 
following relations are satisfied when YTs=i ^ s = (mod/c): 

Se( B I . . . ^.S^n) I Vac)('=) = gfree ( (0, 0) | • • • V'.„K) | (0, 0) ) , (B.35) 

and 

?ie(5|V'e.(^l) ... V'e„(^n)|vac)(^)=gfree WM^e^^,) . . . Ve„ (^n) KMIT • (B.36) 



Here we use the notation | vac)*^'^^ = | (0, 0) ) ® | (0, 0) ) for the closed string ground state of 
the minimal parafermionic model. 

It is essential to point out that correlators ()B.35|) and ()B.36|) are not single-valued func- 
tions of their variables. To resolve the phase ambiguity in ()B.35|) we require that, (a) 



12 



In the case under consideration the modular matrix is given by [10,29], 

(j,M) 2l^^*^'" . ^7r(2j + l)(2J+i; 
b ). ■/ = — , sm ' 



(^■'"^ ^k{k + 2) V fc + 2 



^^Equation HB.33|I requires special consideration for even k and j — j (see Ref. [15] for details) 
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when the arguments are real, Vg = r^, then ()B.35|) is a real function within the domain 
1/T > Ti > . . . > r„ > 0, and (b) the short distance expansions of ()B.35|) must be consistent 
with the OPE's ()2.4|) provided ^'n(T) = V'nl''" ~ ^x)\x=o- Making use of the relation between 
parafermions and SUk{2) WZW currents [5], it is possible to show that both conditions can 
be indeed satisfied. Similarly, one can resolve the phase ambiguity of ()B.36j) . As a result, the 
boundary values of the correlation functions ()B.35j) and ()B.36|) at real Vs = Vg = Tg coincide 
within the domain 1/T > ti > . . . > Tn > 0. 

B.3 Baxter's operators 

In the boundary state formalism, thermal expectation values flA.14|) are expressed in terms 
of correlators of the form: 



free( B 




vac ) 


free( 


B 


vac ) 



with |i?)fi.ee defined in flB.l|) . Our previous analysis implies, therefore, that the integrand 
in ()A.26|) can be written as an expectation value over the highest vectors \ P,{J, M) ) = 
\P) ®\{J,M)) with J = M = 0: 

((0,0),P|U,,(ri)---U,„(T:„)|P,(0,0)) . 
Here the parameter V is related to the zero-mode momentum P: 

0P 

V = 27iiT , (B.39) 



and 



U±(r) = V± e^' ^"^(r) . (B.40) 



From the mathematical point of view, U-|-(f ) are intertwining operators (chiral vertex oper- 
ators) of the conformal dimension d ()A.4|) . They act in the chiral Hilbert space 

Hfree = y" ©(J,M) Vp,(J,Af) with Vp, (J,Af ) = J^p (g) yV( J,M) (B.41) 

in accordance with the formula 

U±(r) : Vp,(j,M) ^ Vpi^_(^j^^ii) . (B.42) 

Equation ()B.38|) . in turn, implies that the power series A±{s,p) ()A.24|) . ()A.30|) can be thought 
of as expectation values over the state | P, (0, 0) ) with P = 2\/kp/(3 of the operators 



uj e vk 



Tr 



X (B.43) 
T2MH ^±2.i^^UoH rexpj/i^'^^dr K±(r)} 
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Here K±(r) = a^U+(r) + a±U^( r) and s = /i^ 

Some explanations of formula ()B.43|) are in order at this point. First, it contains the 
ordered exponential (the symbol T denotes the path ordering) which is defined in terms of 
the power series in /i as follows, 

Texp^fiJ dr K±(r) I = 1 + X (B.44) 

/ dri . . . dr„ K±(ri) K±{t,) . . . K±(r„) . 

Jl/T>Tl>T2>...>Tn>0 

In fact, only even powers of n survive in the definition ()B.43|) . so A± are formal series in s oc yU^ 
whose expansion coefficients are operators acting in the chiral Hilbert space ()B.4H1 . Second, 
the factor u^"^^^ appearing in ()B.43jl was undefined yet. As follows from the commutation 
relations ()A.19|1 . only integer powers of the operator u"^^ are involved in ()B.43jl . The action 
of u"^^ in the chiral Hilbert space ()B.41|) is defined by the formula 

UJ^^ Vp, (j,Af ) = u;^^' Vp, ^j,M) . (B.45) 

Notice that, because of the isomorphism >V(j,_j) — W(fe_jfc_j-), the operator M itself is 
ill-defined. 

In spite of the fact that the factor u"^^^^ does not contribute to the matrix element 
(P, (0,0) I A± I P, (0,0) ), it plays an important role in the definition of A±. Due to this 
phase factor the operators A± commutes with the operator 

li = 2^ ((50)^ + ^-2) . (B.46) 

A proof of the commutation relation 

[Ii, A±(s)] = 0, (B.47) 
is based on the quasiperiodic properties of the intertwining operators if^^{u): 

^A±(m + 1/T) = tl^^iu) u'^^^ . (B.48) 

It follows along the line of Appendix C of Ref . [26] . Notice that the operators Ii + Ii and 
Ii — Ii coincide, respectively, with Hamiltonian and the momentum operator in the closed 
string channel. Hence, commutation relation ()B.47|) manifests the invariance of A±{s) with 
respect to translations along the r-direction. 

It is worth to note that the operators A-|-(s) ()B.43|) act invariantly in Vp,(j,m)- Further- 
more, since each space Vp,(j,m) naturally sphts into the direct sum of finite dimensional "level 
subspaces" 

Vp,(.,M) = ©r=o V^Im) ; ^ ^PUM) = ^ ^PUM) ' (B-49) 

and the grading operator L essentially coincides with Ii, then A±{s) act invariantly in all 
the level subspaces Vpjjj^^-). In particular, the highest vectors | P, {J,M) ) are eigenvectors 
for A±{s) and 

0P 

A±(s)|P,(0,0)) = A±(s,p) |P,(0,0)) with P=f7=. (B.50) 
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Following the route of the work [26] it is possible to prove that A-|-(s) commute among 
themselves for any values of s, 

[A±(s) , A±(s') ] = [A+(s) , A_(s') ] = , (B.51) 

and satisfy the operator-valued identity 

e'^^' A+iqs) A4q-h) - 00^^ A+iq-'s) A4qs) = e'^^' - 00^^ , (B.52) 

with q = e"^"^. Again, formulae ()B.5Hl and ()B.52|1 should be understood as an infinite set of 
formal operator-valued relations for the expansion coefficients of A± without any reference 
to convergence issues. The quantum Wronskian relation ()A.38|) follows immediately from 

(|R50|) and (lR52|) . 

A key idea of the proof of ()B.51|) and ()B.52|) is based on the observation that the formal 
integrals 

1 f^/^ 1 r^/^ 

xq = / du U-iu) , Xi = / du U+(u) (B.53) 

q-q'^ Jo q-q"^ Jo 

satisfy the Serre relations for the quantum affine algebra f/g(sZ(2)) and can be identified with 
the generators of its Borel subalgebra. With this observation the construction elaborated 
in [26] can be applied without any changes. Finally, we note that the operators A-t, acting 
in the CFT chiral Hilbert space, play a role similar to that of the Q-operators of Baxter's 
lattice theory [69]. 

C Appendix : Partition function of the BP sinh- 
Gordon model 

If the parameter (3'^ = —6^ is negative the Schrodinger equation ()A.50|) takes the form ()2.32|) . 
It is quite natural to expect that this differential equation within the parametric domain 
6^ > is somehow related to the BP sinh-Gordon model (j2.14j) . The precise relation is 
proposed in Eq. (|2.36|) . Unfortunately, a rigorous derivation of this formula does not currently 
exist. Here we briefly argue in favour of ()2.36|) based on the high- and low-temperature 
expansions of the partition function Z^l^{h). 

C.l High-temperature expansion 
C.1.1 The BP Liouville model 

First of all, let us discuss the leading high-temperature behavior of ()2.29j) based on properties 
of the Hamiltonian ()2.14j) . Qualitatively, one may expect that if ^eh is not too small, the 

limit T — cxD of (j2.29j) is controlled by either the boundary operator V+ = ^^e^ ^ or 

V_ = ^'_e in (|2.14p depending on the sign ^eh, with some crossover at small 3?e/i. 

More precisely, let us assume that 3?e/i < 0. Then we can treat the term in the 

Hamiltonian ()2.14|) as a perturbation, and in leading approximation one has 

zHMt^oo ^ Z;;\h) {^eh<0), (C.l) 
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where zl^'^\h) is the partition function corresponding to the Hamiltonian of the "BP Liou- 
ville" modeh 

h\^^ = f/free + h^B-l^^+ eTS*^ . (C.2) 



An important feature of the QFT ()C.2|) is the conformal invariance. To be precise, the 
theory can be made scale and conformally invariant by an appropriate redefinition of the RG 
transformation, namely by supplementing it by a formal field redefinition $ ^ $ + ^ ^^rg, 
where the RG "time" ~ ^og{l/T) and Q = b + b^^. In the nomenclature of string theory 
this corresponds to introducing a linear dilaton which modifies slightly the stress-energy 
tensor of the model: 

2^(2) _ _(^d^f^^d2^^W2, (C.3) 



The first two terms in each equation ()G.3|) constitute the corresponding chiral components 
of the modified stress-energy tensor of the Gaussian theory, while W2 and W2 are chiral 
components of the stress-energy tensor of the minimal parafermionic model (see the OPE 

()B.15|) ). The scaling dimension of the boundary operator ^'_|_e^*^ with respect to the mod- 
ified stress-energy tensor ()G.3|) equals one and, hence, the temperature and /i-dependences 
of the partition function zl^\h) readily follows from the dimensional analysis: 

Z^\h) = (27rT)'^ ^-i^ G'('=)(P I b) , (C.4) 

with 

P = i^. (C.5) 

Here G^''^ is some function of the dimensionless parameters P, b and k. 

In a view of the standard heuristic arguments of Ref. [72], it is expected that G^''\P \ b) 
is an analytical function of the complex variable P in the lower half plane QmP < 0. 
Furthermore, following [72] in the "perturbative" calculation of the partition function 
one can first integrate out the constant mode of the Bose field This integration produces 
simple poles in the variable P at the points P = in by/k (n = 0, I ■ ■ and the corresponding 
residues 

G^:^ = Resp^,,^ „ [ GW(P I 6) ] (n = 0, 1, 2 . . . ) (C.6) 
are expressed through the integrals of the k x n-points Matsubara correlation functions 

27ri 

r dn... dr,„ (( vr^^Vi) • • • vf^'''\n,,) ))'o . 

Jo Jo 

Here gfree is the boundary degeneracy of the free CBC in the minimal parafermionic model 
(see Eq. (II .21) ). and the prime means that the constant mode contribution is excluded from 
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the thermal averaging (( ■ ■ ■ ))o ()A.14|) . The free- field correlators in ()C.7|) have an especially 
simple form for n = 1. In this case the integral ()C.7|1 can be brought to the form of the 
Selberg integral [63,64]. This yields the explicit formula 



iD gfrcc 



2vri kz 

gP gfrec -j- 

27ri 







/ dvi... dVk 


Jo 


Jo 




271 




[v^r(i-i±^)J 



d.;. n [2 sin (^) 



2(1+6-') 
k 



k 



(C. 



Another highly nontrivial property of the BP Liouville partition function can be guessed 
by virtue of the form of stress-energy tensor ()C.3|) . Namely, it is invariant with respect 
to the transformation h h^^ . In the case of the conventional Liouville model a similar 
phenomenon manifests a remarkable non-perturbative duality of the theory. If one admits 
that the same symmetry occurs in the BP Liouville model, then, the function G^^'' (P | h) 
should also possess an additional series of "dual" poles at P = ih~^\/k n. 

It should be mentioned that the model ()C.2|) has been already studied in the works [70] 
and [71] for the cases k = 1 and k = 2, respectively. In both cases the function G^'''^ ()C.4|) 
looks as follows 



G^''\P\b) 



gPgfrcc 

27riP 



r 1 + 



ibP 

Vk 



r 1 + 



bVk^ 



27T b- 



yfcr(i-i±^; 



iP^/k 



(C.9) 



Notice that ()C.9|1 satisfies the all above-mentioned conditions even thought k > 2. 



C.1.2 General structure of the high-temperature expansion 

Above we have discussed the leading high-temperature behavior of ^^shg in the case 'Steh <Q. 

Obviously for 3?e /i > 0, z!^\h) in asymptotic formula ()C.1|) should be replaced by zl^\—h). 

Since zl^\±h) vanishes in the limit T ^ 00 if /i is taken in the "wrong" half plane (note 

the factor T in ()(I4jl ). this in turn suggests that the overall T ^ 00 asymptotics of the 
partition function is correctly expressed by the sum 

4ilT-.oo-4f^W+4fV/^)- (CIO) 

What can be said about corrections to the leading asymptotic? Again we assume that 
^eh < and consider the perturbative effect of the term fi V- in (I2.14j) to the partition 
function ()C.4|) . In the unperturbed BP Liouville model theory the parameter /i is a dimen- 
sionless constant. Let us eliminate n from the Hamiltonian ()C.2|) by shifting the field 
Then the coupling fi in front of the in Eq. ()2.14|1 is replaced by /i^. Since the anomalous 
dimension of the boundary operator V_ with respect to the BP Liouville stress-energy ten- 

2bQ 

SOT ()(7.3|1 is given by 1 - then ~ £'* , where E^: is the physical energy scale of the 
BP sinh-Gordon model. Hence we deduce that the perturbative corrections to the leading 
asymptotics ()C.1|) should be in a form of power series expansion of the dimensionless param- 
eter K~ , with K = E^/T. The case 3?e /i > can be analyzed similarly and one comes to 
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the same conclusion about the form of perturbative corrections to ()C.10|) . As was mentioned 
above, it is quite natural to expect that the leading asymptotics ()C.10|) is invariant with 
respect to the duality transformation b —>■ b^^. If one assumes that the BP sinh-Gordon 
partition function possesses this property of self-duality as well, then the high-temperature 

^bshg 



(k) 

expansion of ^bsL should be of the form 



= Zi!\h) MiK, h) + Zl^^\-h) MiK, -h) . (C.ll) 



bshg 

2bQ 2Q 

Here M{k, h) is a formal double power series in integer powers of k ^ and k bk ; 



M{K,h) ~ l + ^M„,„(/l) f^T-i^b+mb'^) _ 



C.1.3 Small k expansion of the spectral determinant 

In order to justify the proposed formula for ^^shg (ESSJ) we should observe the structure 
(ICIllj) in the small-K expansion of the spectral determinant D{k,,C,) ()2.42j) . When k, goes 
to zero and u ^ 1 the potential in the Schrodinger problem ()2.H2|) can be approximated by 

„ ku 

K e'>Q . To understand the quality of this approximation one can make a change the variable 
()A.57|) with d = 1 — bQ/k and bring the differential equation (j2.32j) to the form ()A.58|) . 
Therefore for u ^ 1 

. (C.13) 

V ^/k 

Here and below, in order to write formulae in the most instructive form, we use the notation 
P = ^ ^ rather than ^ (see relations ()C.5|) and ()2.38|) ). The overall normalization in ()C.13|) 
is chosen to ensure the asymptotic condition ()2.34j) . Within the domain, 

-i ^og{j,)«u«flog{^), (C.14) 

the potential in ()A.58|1 develops a wide plateau and the solution 0+ ()(113j) becomes a 
combination of the two plane waves. 



0+H = V^^r -1— — - A{K,-P\b)e^Q + 



27r/c V V k 

with 

2Qb 



27Tk V y/kJ V k 

^rfi^) f^)^^A(.,P|6)e-^^ (C.15) 



A{K,P\b) = 1 + 0{k—) . (C.16) 
When K is small but finite, corrections to ()C16|) can be obtained using the perturbation 



theory. In view of equation ()A.59|) . these corrections have the form of a power series in k ? : 

00 

A(k,P|6) ~ 1 + ^ A„(P|6) . (C.17) 



2Qbn 

n=l 
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The formal power series A{k, P \ b) amounts to A^{s,p) ()A.55|) taken at = —6^ < with 
the variable s related to k, as in ()A.56|1 and p = i^^- We have no reason to expect the 
convergence of the power series ()C.17|) for 6^ > 0, i.e., the series should be understood as a 
formal asymptotic expansion. 

To find the form of the solution B_ in the domain ()C.14|) . we note that the transformation 
u —>■ —u and h leaves the Schrodinger equation ()2.H2j) unchanged while interchanging 

the solutions 0+ and B_ ()2.H4j) .^^ Hence, 0_ can be obtained by means of this transfor- 
mation from ()C.15|) within the domain ()C.14|) . Now the Wronskian ()2.42p can be calculated 
at any point from ()C.14|) . where both Q± are combinations ofplane waves. It yields the 
following form of the small-K expansion of the spectral determinant: 

D{k, ^ RiP) A{k, P I h) A{k, P I h-^) + R{-P) A{k, -P \ b) A{k, -P \ b'^) , 

with 



iP fibP\ f iP \ ff^Q 



andP = ^e 

The high-temperature expansion (jC.18|l is in agreement with (IC.llll . Furthermore, it 
implies that formula ()(y.9|l holds true for any k > 1 . 



C.2 Low-temperature expansion 

Here we elucidate the low-temperature behavior of ^bshg(^) using the boundary state for- 
malism described in Appendix B. 

In the case of the BP sinh-Gordon model the boundary state | P)bshg is some particular 
vector in the closed string space of state 

|5)bshgG J ®iJ,M)^PUM)®^P,iJ,M) ^ (C.19) 

where we use the same notations as in formula ()B.4H) . Unlike the conformal boundary 
states discussed in Appendix B, | B )hshg essentially depends on the temperature via the 
dimensionless parameter k = E^^/T. In the presence of the external field h ()2.14|) . the 
partition function ^bshg(^) can be expressed in terms of the vacuum overlap of the boundary 
state ()(y.l9|l analytically continued to pure imaginary values of the variable P ()(I5|) : 

^bshg(/i) = bshg( P I ■ { I P, (0, 0) ) ® |P, (0,0)) } 

The motivation behind ()C.20|) can be found in Ref. [1] where a similar relation was exploited 
in the context of the circular brane model. 

"'^^For this reason the partition function l|2.36|) is invariant with respect to the duahty transformation 



(C.20) 
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The simplest idea about the infrared fixed point of the BP sinh-Gordon boundary fiow 
is that it corresponds to the fixed CBC for both bosonic and parafermionic sectors of the 
theory. This scenario imphes that 

|5)bshg ^e"^ lE)!^*^ ® |5(o,,„)) as K^oo. (C.21) 

Here -Ebshg oc E^. is the ground state energy, | B(Q^rn) ) is the boundary state corresponding to 
the fixed CBC of the minimal parafermionic model ()B.31jl .^^ and 

/oo 
n=l 

is the boundary state associated with the Dirichlet CBC [67,68]. The assumption (jC.2H) 
leads, in turn, to the leading low-temperature asymptotics for the BP sinh-Gordon partition 
function: 

logZbshg(/^) = + log (gDgfixed) + 0(T) . (C.23) 

To make more definite predictions about the low-temperature expansion of ^bshg(^); one 
should take into account the integrability of the model. In a forthcoming paper [73] we intend 
to present a comprehensive discussion of a relevant set of local Integral of Motions (IM's) 
in the context of more general QFT model with boundary interaction. Skipping details of 
the computational complexity we mention here that the boundary state ()(I19jl is expected 
to satisfy an infinite set of conditions 

( I21-1 - I21-1 ) I B )bshg = ( / = 1, 2 . . . ) (C.24) 

for some operators 

l2i-i{x) = —T,i{r- ix) , hi-i{x) = / — T2i{t + ix) , (C.25) 
Jo 27r Jq I-k 

where T21 and T21 are chiral currents of the Lorentz spin 21 and (—2/), respectively. As 
has been explained in Appendix C, the local chiral currents of the parafermionic minimal 
model constitute Wfc algebra. Therefore the local densities T21 are, in fact, appropriately 
regularized polynomial in dcf) and Wjt currents as well as their derivatives. Similarly, T21 are 
constructed using d(f) and currents. Notice that any operators in the form ()C.25|) are 
conserving charges in the sense that 

^l2/-i(x) = ^i2/-i(x)=0, (C.26) 

and the meaning of ()C.24|) is that the boundary neither emits nor absorbs any amount of 
the combined charges "^21-1 — "^21-1 [55]. For this reason I21-1 and I21-1 are referred to as local 
IM's. 



^^Recall, in order to fix the value of integer m = 0, — 1 in (|C.21() . one needs to specify the 

normahzation of the boundary fields *± unambiguously (see Section |OJ. 



49 



The local IM's Ii and I3 are known in the closed form. Namely, Ii is given by ()B.46|) and 
an explicit bosonized form of the first nontrivial local density T4 is presented in Ref. [74]. 
There are strong indications that the commutativity conditions 

[hi-i,li] = [hi-i,h] = (C.27) 

(a) fix all the operators I21-1 for / > 2 uniquely up to the normalization, and (b) such defined 
operators I21-1 mutually commute: 

[l2,„i, Wi] =0 . (C.28) 

In order to specify l2«~i unambiguously, it is convenient to normalize the corresponding local 
density in such a way that the monomial (50)^^ appears in T21 with the coefficient equals 
one, i.e. 



I 



21-1 



We can now take advantage of the general prediction from the Introduction of Ref. [31] 
about the form of the low-temperature expansion in the case of an RG boundary flow which 
possesses an infinite set of mutually commuting local IM's, and terminated at the "triv- 
ial" infrared fixed point (the one which in Cardy's classification [13] corresponds to the 
identity primary state). In the case under consideration the prediction implies that the 
low-temperature expansion of logZbshg(^) looks as follows 

logZbshg(/i) ^-fo Y + log (gi?gfixed) -2^fi hi-i[i , (C.30) 

1=1 

where /; are some coefficients and l2i-i{P) are vacuum eigenvalues of the operators ^21-1'- 

I21-1 I P, (0, 0) ) = {27iTf-^ l2i-i{P) I P, (0, 0) ) . (C.31) 

Let us emphasize that the l2i-i{P) are unique polynomials of order / of the variable P^, and, 
as follows from the normalization condition ()C.29|1 

l2i-i{P)=(^) +... , (C.32) 



2 

where the omitted terms contain lower powers of P^. In particular, in view of the explicit 
formulae for Ii ()B.46|) and I3 [74] it is possible to show that 

- Q^-w^)^ ^^-^^^ 

' \2J 4{2 + 3k)\2J 64 (2 + A;) (2 + 3A;) ' 

with Q = b + b~^. 

Applying the standard WKB iterational scheme [33] to the spectral determinant (j2.42p . 
one can reproduce exactly the structure of the low-temperature expansion ()C.30|) as well as 
the eigenvalues ()C33|) . En passant, the WKB calculation provides an explicit form of the 
the expansion coefficients: 

r(/-i) r{^b)T{^-^b-^) fVkY 
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D Appendix : Low-temperature expansion of ' 

At low temperature the regularized partition function ()2.53|) can be studied using the WKB 
expansion. The leading WKB approximation yields an explicit expression for the regularized 
ground state energy ()2.56|) : 



E. 



(fc) 



Tim — 



2 (1 - i) + E T7?i^ [v^ (- - i) - ^ (-^ - 1) + ^ (I) - (I) J -n-(^) . 

m=2 V 



(D.l) 



It is useful to bear in mind that ()D.1|) is applicable for < < | only. It can be analytically 
continued within the domain | < ^ < vr by means of the relation: 



1=1 



I sin2 (^(^) (1 _ cos(^)) log [sin (§)] + (l + cos(^)) log [cos (§)] 



(D.2) 



The following formula for the regularized ground state energy holds within < ^ < vr 

EfV^* = C^f -(7i? + 21og2 + ^(|)-logA;) cos(e) + (D.3) 



where 



Vtt 2 

k <=os 



1,11 



- - 1-2 - 



COS 



\e)) 



h CO^'W Ei 3^2^ + 1, 1,1; 2,11 COS^(0)) 



1=1 



(fc) 



fc-1 _ l^ 

>P ^ U 2^ 

1=1 ^\k) 



gi^2 



(D.4) 



The systematic WKB expansion of ()2.53p leads to an asymptotic series of the form 



logZ, 



(k) 



E, 



(k) 



T 



+ l0g(gfixed) + E^/'^(^) 



1=1 



T \ 2^-1 



(D.5) 



where the coefficient Fi''^ reads explicitly as 



k-1 



12{k + 2) 
24r(i-l) L ^ 



2^(l4 + i;l + ilsin^W) 

2 A;- 1 



(D.6) 



12A; 



cos(e)2Fi(l + i,l;| cos2(^^)) 
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